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I^d  ABSTRACT  (MMfmum  iOO  wordsj 

The  purpose  of  this  experiment  was  to  perfoim  a  controlled  field  experiment,  involving  the  injection  of 

several  aromatic  hydrocarbons  and  a  nonreactive  tracer  into  an  uncontaminated  aquifer.  By  monitoring 
the  plume  development  of  these  solutes,  and  by  measuring  a  number  of  physical  and  chemical 
characteristics  of  the  aquifer,  this  study  was  deigned  to  provide  data  on  those  properties  which 
significantly  control  the  propagation  of  dissolved  contaminants  in  groundwater  systems.  A  secondary 
objective  was  to  measure  the  in  situ  degradation  rates  of  the  selected  organic  conqiounds.  The 
dis^ipearance  and  transformation  of  the  organic  solutes  during  this  experiment  demonstrated  that  natural 
d^radation  processes  were  able  to  effectively  reduce  these  levels  of  dissolved  organic  contaminants  in  a 
reasonable  time  frame.  This  represents  the  most  important  result  of  this  field  study,  since  active 
remediation  would  not  be  needed  in  situations  where  natural  degradation  rates  were  sufficiem  to  reduce 
contaminant  concentrations  to  safe  levels.  The  monetary  and  environmental  cost  savings  of  allowing 
natural  biological  restoration  of  the  residual  contaminants  are  potentially  enormous.  It  appears  that  tiiis 
study  is  the  first  field  oqieriment  to  prove  conclusively  that  hydrocarbon  solute  losses  were  due  to 
chemical  degradation  rather  than  physical  losses. 
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OV  ANSi  Sto 


PREFACE 


This  report  was  prepared  by  the  Armstrong  Laboratory  Environics  Directorate  (AL/EQC), 

139  Bunes  Drive,  ^ite  2,  Tyndall  Air  Force  Base,  Florida  32403-5323,  and  the  Tennessee 
Valley  Authority  (TV A)  &igineering  Laboratory,  Norris,  Tennessee  37828,  in  conjunction 
with  the  Electri^  Power  Research  Institute  (EPRI).  The  reported  work  was  cofiinded  by  the 
Air  Force  aixl  EPRI  through  an  interdepanmental  procurement  action  with  TVA. 

This  final  report  describes  the  second  MacroDispersion  Experiment  (MADE-2)  conducted  at 
Columbus  AFB,  MI;  the  experimental  methodologies  that  were  used  to  analyze  site 
characteristics,  inject  chemical  tracers  and  solutes  into  the  aquifer,  and  monitor  the  resulting 
plume;  the  data  analysis  techniques  that  were  used;  and  an  explanation  of  the  significance  of 
these  experimental  observations.  EPRI  TR-101998  provides  a  more  general  description  of 
the  experiment. 

The  work  was  performed  between  October  1989  and  September  1992.  The  AL/EQC  project 
manager  was  Dr.  Thomas  B.  Stauffer. 
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EXECUTIVE  SUMMARY 


A.  OBJECTIVE 

The  purpose  of  this  experimeot  was  to  perform  a  controlled  field  experiment, 
involving  the  injection  of  several  aromatic  hydrocarbons  and  a  nonreactive  tracer  into  an 
uncontaminated  aquifer.  By  monitoring  the  plume  development  of  these  solutes,  and  by 
measuring  a  number  of  physical  and  chemical  characteristics  of  tt^  aquifer,  this  study  was 
designed  to  provide  data  on  those  properties  which  significantly  control  the  propagation  of 
dissolved  contaminants  in  groundwater  systems.  Data  of  this  type  are  quite  rare,  and  are 
extremely  useful  for  validating  numerical  fate  and  transport  models  which  include  complex 
processes,  such  as  solute  advection-dispersion,  geochemical  attenuation,  and  microbial 
degradation.  A  secondary  objective  was  to  measure  the  in  situ  degradation  rates  of  the 
selected  organic  compounds,  in  order  to  investigate  the  possibility  of  utilizing  natural 
attenuation  processes  as  an  alternative  to  active  ren^iation  techniques. 

B.  BACKGROUND 

Groundwater  contamination  is  an  important  environmental  concern  for  both  military 
and  civilian  industrial  operations.  The  organic  constituents  of  fuels  such  as  JP-4,  gasoline, 
and  diesel  (e.g.,  benzene,  toluene,  nafrtithalene,  xylene,  etc.)  and  cleaning  solvents  (e.g., 
trichloroethylene)  have  resulted  in  serious  contamination  problems  at  many  DOD  facilities. 
Many  of  thke  chemical  components  of  fuels  and  solvents  are  regulated  by  EPA  as  hazardous 
substances,  and  have  stringent  concentration  limits  in  drinking  water  (Geraghty  and  Miller, 
Inc.,  1991).  Their  presence  in  groundwater  at  concentrations  exceeding  the  ^A’s  maximum 
contaminant  levels  generally  requires  the  use  of  costly  water  purification  methods.  The  lack 
of  an  adequate  quantitative  understanding  of  the  processes  that  govern  the  environmental  fate 
of  these  contaminants  has  constrained  the  development  of  cost-effective  groundwater 
remediation  techniques. 

A  great  number  of  research  efforts  have  studied  various  physical,  chemical  and 
biological  mechanians  that  influence  the  transport  of  pollutants  in  the  environment.  Most 
have  been  performed  within  laboratories,  but  some  relatively  recent  experiments  have  been 
conducted  at  fkld  sites  to  determine  the  effects  of  natural  conditions  and  heterogeneities  on 
these  processes  (Boggs  et  al.,  1992;  Harvey  et  al.,  1993;  Sudicky,  1986).  Additionally,  the 
accidratal  discharge  of  contaminants  into  groundwater  systems  h^  allowed  some  limit^ 
studies  of  die  in  situ  transport  of  pollutants.  However,  the  lack  of  specific  information  about 
the  initial  conditimn  of  contaminant  mass  and  location  hinders  the  ability  to  draw  meaningful 
conclusions  from  such  research. 

Previously,  a  few  field  experiments  have  reported  the  attenuation  of  test  solutes,  and 
claimed  diat  the  results  were  possibly  due  to  microbial  degradation  processes  (e.g.,  Sutton 
and  Barker,  1985).  While  laboratory  studies  have  proven  that  soil  organisms  are  capable  of 
utilizing  a  great  variety  of  chemical  compounds,  similar  work  showing  natural  degradation  in 
the  field  has  not  previously  been  possible,  for  a  number  of  reasons.  Although  it  is  generally 
accqited  that  the  in  situ  biodegra^tion  of  chemical  contaminants  does  occur,  conclusive 
experimental  evidence  is  lacking  in  the  current  literature. 
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c.  score 


tepoct  describes  the  research  conducted  at  a  field  site  at  Coiumbus  Ak  Force 
Base,  Miaiss^i.  Section  I  provides  an  introductkm  to  the  objectives  of  this  study  and 
pertioeiit  background  irdbnnation.  Section  n  describes  the  physical  and  geochemi^ 
dMuicteristics  of  die  test  site.  Section  in  outlines  the  meth^  that  were  used  to  conduct  the 
field  eiqieriinent.  Section  IV  discusses  the  results  of  the  data  analysis  and  some  general 
observatHHis.  Section  V  presents  a  more  detailed  madiematical  overview  of  the  ^tial 
nKMneitts  analysis  of  the  solute  plumes.  Section  VI  shows  how  the  degradation  rates  were 
calculated  for  die  aromatic  hydrocarbons.  Section  vn  reviews  the  dispersivity  computations, 
and  section  VIH  presents  the  researchers  concluskms  and  recommendations. 

D.  METHODOLOGY 

To  effectively  monitor  the  development  of  the  plume  of  hydrocarbons  and  die  tritium 
tracer,  a  well  field  consisting  of  328  multilevel  and  56  BarCad*  samplers  was  available  to  the 
researchers.  Because  of  time  and  cost  constraints,  only  five  complete,  three-dimensional 
sanqile  sets  were  taken  at  about  lOO-day  intervals  during  the  study.  Additional,  smaller 
sampUng  events  were  used  to  observe  phime  activity  on  a  shorter  time  scale.  Since  it  is 
neiter  jdiysically  nor  economically  possible  to  sanqile  an  entire  aquifer,  a  three-dimensional 
network  of  wdls  must  be  used  to  determine  the  location  and  concentration  of  groundwater 
contaminants.  To  obtaui  a  conqilete  "picture"  of  this  plume,  geostatistkal  techniques  (e.g., 
kriging)  were  enqdc^ed  to  mathematically  fill  in  die  spatial  gaps  between  sampllrq^  wells. 
nUs  techni^  is  the  most  sc^ihisticated  and  accurate  medxxl  available  for  use  in  generating 
stochaittic  niodels  and  conqaiter  grafdiics  of  subairfM:e  contaminant  transport. 

Aldiough  die  measuranent  of  degradatkm  rates  was  a  principal  objective  of  diis  smdy, 
no  attempts  were  made  to  classify  or  quantify  the  microbial  populations.  While  a  detailed 
analysis  of  mzymatic  activity  before,  during  and  after  the  nrovement  of  the  plume  through  a 
given  area  would  greatly  support  observations  of  biodegradation,  such  work  was  beyond  the 
scope  of  this  study  and  is  left  for  ftmire  research  efforts.  In  order  to  distinguish  solute 
degradation  from  losses  due  to  sorption,  dilutkm  or  evaporation,  and  to  permit  the  detection 
of  dqpulation  products,  ’^-labeled  p-xylene  was  mctuM  in  the  injection  solution. 

Frevious  work  at  this  site  characterized  many  of  the  important  physical  pn^ierties  of 
this  aquifer  (Boggs  et  al.,  1992).  In  this  study,  soil  samples  were  taken  at  various  locations 
and  d^pdis,  then  analyzed  for  particle  sizes,  organic  carbon  content,  surface  area,  and 
partilkm  coefficient  Grc^water  sanqiles  also  tested  for  a  number  of  elemental 
components  and  chemical  characteristics.  Finally,  a  set  of  piezometers  was  used  to  measure 
tenqpoiil  floctuatic^  in  the  hydraulic  gradiem  of  the  aquifer,  which  were  due  to  a 
conAinatkm  of  both  daily  rain  events  and  evapiHranspiration,  and  seasonal  factors. 

E.  TEST  DESCRIPTION 

The  natural  gradient  field  experiment  was  omducted  to  investigate  the  tran^rt  and 
degradation  of  four  dissolved  aromatic  organic  conqxxinds  (benzene,  naphthalene,  p-xylene, 
and  o-dichlowbeiMHme)  and  one  nonretained  tracer  in  an  aquifer.  The  ^udy  was  done  at 
CdunduM  AFB,  Minissippi,  in  a  shallow,  unconftned  alluvial  aquifer  having  conskterable 
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heterogeneity  in  physical  and  chemical  properties.  A  two-day  pulse  of  9.7  cubic  meters  of 
tritiated  water  and  the  hydrocarbons  in  a  dilute  aqueous  solution  was  injected  into  the 
aquifer.  The  initial  concentrations  of  the  organic  test  compounds  in  solution  ranged  from 
about  7  to  70  mg/L,  which  was  prepared  and  stored  on-site  using  ambient  groundwater  from 
an  upgradient  local  well.  The  injection  was  made  into  five  wells,  spaced  1  meter  apart  and 
aligiied  transverse  to  the  direction  of  groundwater  flow.  Each  injection  well  was  screened 
over  a  0.6  meter  interval,  which  was  placed  within  the  saturated  zone  of  the  aquifer. 

The  solute  concentration  distributions  were  monitored  at  one  to  three  month  intervals 
over  a  15-month  period,  using  a  three-dimensional  network  consisting  of  over  6000  sampling 
points.  Selected  points  within  the  network  were  sampled  more  frequently  to  produce  time- 
series  sets  of  solute  concentrations.  Spatial  moments  calculated  from  the  three-dimensional 
netwoiic  concentration  data  and  from  the  time-series  concentration  data,  were  used  to 
determine  the  fate  and  transport  behavior  of  each  solute. 

Stable  organic  compounds  in  groundwater  samples  were  analyzed  by  gas 
chromatography  (GC).  Deuterated  toluene  and  4-bromofluorobenzene  internal  standards 
were  added  to  20  mL  water  samples,  which  were  then  extracted  with  2  mL  of  n-pentane. 
These  extracts  were  then  aiudyz^  for  benzene,  p-xylene,  mq>hthalene  and  o-dichlorobenzene, 
also  by  GC.  Sensitivity  of  the  method  was  50  |ig/L  for  benzene  and  4  ^g/L  for  the  other 
conqxMnds.  All  sanq)les  were  extracted  and  analyzed  within  21  days  of  collection.  Quality 
coittrol  for  the  organic  analyses  was  ensured  by  using  several  techniques,  including  the  use  of 
sample  duplicates,  certified  standard  control  solutions,  internal  standards  for  all  analyses,  and 
independent  verification  of  our  analytical  results  by  another  laboratory. 

Tritium  and  in  samples  from  the  wells  were  analyzed  by  liquid  scintillation 
counting  in  dual  istMcpe  mode,  by  the  Water  Resources  Research  Center  at  Mississippi  State 
University.  Samples  were  counted  for  20  mimites  or  to  a  1  percent  error  at  the  95  percent 
probability  level,  whichever  was  attained  first.  Background  levels  of  tritium  and  in 
ambient  groundwater  at  the  test  site  dictated  the  analytical  sensitivity  for  these  measurements, 
and  were  spproximately  2  and  3  pCi/mL,  respectively.  The  statistical  errors  associated  with 
20-tninute  measurements  at  these  levels  are  about  ±20  percent  at  the  95  percent  probability 
level.  A  set  of  tritium  and  '^C  standards  was  counted  every  108  samples  to  check  analytical 
accuracy  and  verify  instrument  performance.  Every  25'*'  sample  was  counted  twice  to 
estimate  meamrement  precision  and  repeatability.  As  an  overall  check  on  the  radiological 
measurements,  duplicate  field  samples  amounting  to  roughly  5  percent  of  the  total  field 
samples  were  analyzed  independently  by  the  TVA  Western  Area  Radiological  Laboratory. 

Dissolved  oxygen  (DO)  measurements  were  performed  at  2-  to  3-month  intervals  at 
selected  locations  during  the  experiment.  DO  data  were  intended  to  demonstrate  whether 
aerobic  conditions  for  possible  degradation  of  the  organic  compounds  were  maintained,  and 
were  collected  from  16  multilevel  samplers  (MLS)  located  along  the  longitudinal  axis  of  the 
tracer  plume,  and  6  piezometers.  Samples  were  collected  in  60  mL  BOD  bottles  after 
purging  ai^roximately  100  mL  from  the  MLS  tubes,  and  measurements  were  made  with  a 
calibrated  DO  probe  irrunediately  after  sample  collection.  The  DO  measurements  at  the  6 
piezometers  were  made  in  conjunction  with  pH,  conductivity,  and  oxidation-reduction 
potential  measurements. 
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F.  RESULTS 


The  tritium  phime  was  characterized  by  extreme  skewness  in  the  longi&Klinal 
threctkm.  The  evolution  of  the  plume  was  consistent  with  the  measured  hydraulic  head  and 
conductive  fields,  and  was  primarily  due  to  the  increasing  mean  hydraulic  conductivity  with 
diriance  frra  the  ii^tion  site.  Measurenwnts  of  hydraulic  conductivity  produced  a  mean 
Darcy’s  Law  groundwater  velocity  of  5  m/yr  near  the  injection  wells,  whereas  local 
velocities  in  the  far  field  were  foiind  to  exceed  400  m/yr.  The  dispersivity  of  the  alluvial 
aquifer  was  estimated  by  fitting  the  tritium  data  to  general  spatial  moments  equations  for 
twcHdinoensional  advective-dispersive  transport  in  a  nominiform  flow  field.  A  longitudinal 
di^)ersivity  of  approximately  10  meters  provided  the  best  representation  of  the  first  and 
secmid  ItH^itudin^  moments  of  the  tritium  plume. 

The  aromatic  hydrocarbons  degraded  significantly  during  the  experiment.  Analysis 
for  die  degndatkm  inducts  of  ^*C  p-xylene  revealed  that  about  80  to  W  percent  of  the 
preseitt  was  associate  with  dissolv^  and  intermediate  products,  indicating  aerobic 
degradation  of  the  p-xylene.  Dissolved  oxygen  in  the  pulse  plume  maintained  aerobic 
conditkms  throughwt  the  experiment,  and  was  always  greater  than  2.6  mg/L.  Degradation 
kinetics  calculated  from  the  whole-field  data  set  were  ai^roximately  first-order  with  the 
following  rate  constams:  benzene,  0.0070  d*';  p-Xylene,  0.0107  d  ';  mqihthalene,  0.0064  d'*; 
OHlichlotobenzene,  0.0046  d  ' .  Similar  reaction  nues  were  also  obtain^  from  a  near-field 
subset  of  the  data,  using  a  model  based  on  the  hydrologic  characteristics  of  the  aquifer.  The 
shiqies  of  die  degradation  rate  curves  were  consistent  with  microbial  degradation  processes. 
Maximum  d^ndatkm  rates  obtained  are  presumed  to  be  characteristic  of  the  microbial 
pt^wlation’s  metab(dtsm.  A  least-squares  inversion  method  of  fitting  time-series  data  to  an 
analytical  one-dimensional  transport  model,  containii^  first-order  decay  and  linear 
equUitnium  sorptkm  ^ms,  yielM  the  followup  mean  first-order  rate  constants:  benzene, 
0.010  (f’;  mqitttludene,  0.013  d*’;  p-xylene,  0.016  d'  and  o-dichlorobenzmie,  O.OQS  d'^ 

Sotptkm  had  some  affect  <«  die  transport  of  die  organic  compounds,  but  it  was  a 
rdadvely  mintnr  process  conqiared  with  the  degradation.  Because  of  the  dmniiiaitt  influence 
oi  the  degradadve  processes,  field-average  retardation  factors  could  not  be  estimated  fnmi  a 
conqMoistm  of  die  mean  displacement  rates  of  the  anmiatic  cmnpound  and  tritium  plumes. 
However,  analysis  of  tune-series  data  at  each  of  17  sarnie  points  located  within  a  distance  of 
tqquoxinudely  20  meters  of  the  injection  point  {uoduced  mean  retardation  factors  for 
bmizrae,  nafdithalene,  p-xylene,  and  o-dichlorobeozene  of  1.20,  1.45,  1.16,  and  1.33, 
respecdvety.  Agremnem  between  these  results  and  retardation  factors  calculated  from 
sorption  coefficients  (ditained  in  laboratory  batch  and  colunm  experimems,  suggests  that 
lahmatofy  retardation  measurements  may  be  reliable  indicators  of  sorptitm  behavior  in  the 
field. 

G.  CONCLUSKWS 

The  ^Hsappearance  and  transformation  of  the  organic  solutes  during  this  experiment 
demonstrated  thin  natural  degradation  processes  were  able  to  efiectively  reduce  these  levels 
of  dtsstdved  osga^  contaminants  in  a  reasonable  time  frame.  This  result  suggests  that,  for 
rimdar  orgi^  conqxNinds,  it  might  be  best  to  reririct  aquifer  remediation  activities  to  the 
contaminant  source  r^km.  By  reducing  this  source,  the  resulting  plume  of  or^inic  solutes 
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could  possibly  maintain  a  steady-state  limit,  given  the  correct  physical,  chemical  and 
biological  conditions.  In  such  an  aquifer,  the  spatial  boundary  of  this  plume  would  be 
determined  primarily  by  the  hydrology  and  geochemistry  of  the  site,  solute  sorption, 
biodegradation  by  indigenous  microbes,  redox  capacity,  and  oxygen  and  nutrient  supply. 

This  observation  represents  the  most  important  result  of  this  field  study,  and  could 
have  in^rtant  inqilications  for  the  restoration  of  aquifers  contaminated  by  organic 
chemicals.  Active  remediation  would  not  be  needed  in  situations  where  natural  degradation 
rates  were  sufficient  to  reduce  contaminant  concentrations  to  safe  levels  before  they  were 
transported  to  regions  where  they  might  be  dangerous  to  human  health  or  the  environment. 
Rem^iation  by  natural  degradation  processes  is  likely  to  be  applicable  to  residual  organic 
contaminants  (i.e.,  NAPL  organics  immobilized  by  capillary  forces  that  cannot  be  effectively 
removed  from  the  aquifer  by  pumping).  The  monetary  and  environmental  cost  savings  of 
allowing  natural  biological  restoration  of  the  residual  contaminants  are  potentially  enormous. 
After  a  review  of  the  current  scientific  literature,  it  appears  that  this  study  is  the  first  field 
experiment  to  prove  conclusively  that  hydrocaibon  solute  losses  were  due  to  chemical 
degradation  rather  than  physical  losses,  such  as  sorption  and  vaporization. 

H.  RECOMMENDATIONS 

Further  research  is  needed  to  build  upon  results  of  this  study,  to  further  support  the 
observation  of  natural  attenuation  and  to  attempt  to  quantify  the  environmental  requirements 
for  this  process.  The  degradation  rates  determined  for  the  four  aromatic  compou^s  in  the 
Columbus  aquifer  are  now  being  used  in  the  design  and  modeling  stages  of  a  new  field  test 
which  will  use  an  emplaced  NAPL  source.  This  test  at  the  Columbus  site  will  be  a 
demonstration  of  the  practicality  of  natural  degradation  (natural  attenuation)  as  a  remediation 
action  for  a  steadily  leaching  source  in  an  aquifer.  This  experiment  more  closely  emulates 
the  effects  of  a  fuel  spill.  Assuming  that  a  steady-state  contaminant  concentration  situation  is 
attained  in  this  test,  natural  degradation  will  be  a  verified  groundwater  contaminant  treatment 
cation. 

It  is  proposed  that  experiments  similar  to  this  one  should  be  done  at  several  other  test 
sites  to  prove  the  general  applicability  of  this  method,  and  to  measure  both  aerobic  and 
anaerobic  degnulation  rates  for  use  in  the  development  of  more  accurate  fate  and  transport 
moctels.  The  technique  of  using  ‘^-labeled  compounds  is  suggested  for  future  field 
experiments,  because  it  distinguishes  solute  degradation  from  solute  losses  by  sorption  and 
eviqKnatkm,  and  permits  the  demonstration  of  more  accurate  mass  balances  throughout  the 
course  of  the  investigation.  If  multiple  test  sites  become  available,  they  could  then  be  used 
to  conduct  NAPL  studies  which  mi^t  further  support  the  natural  attenuation  treatment 
option,  and  might  even  suggest  this  technique  as  die  action  of  first  choice  in  certain 
groundwater  contamination  situations. 
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SECTION  I 
INTRODUCTION 


OBJECTIVE 

The  purpose  of  this  experiment  was  to  perform  a  controlled  field  experiment, 
involving  the  injection  of  several  aromatic  hydrocarbons  and  a  nonreactive  tracer  into  an 
unctmtaminated  aquifer.  By  monitoring  the  plume  development  of  these  solutes,  and  by 
measuring  a  numter  of  physical  and  chemical  characteristics  of  the  aquifer,  this  study  was 
designed  to  provide  data  on  those  properties  which  significantly  control  the  propagation  of 
dissolved  contaminants  in  groundwater  systems.  Data  of  this  type  are  quite  rare,  and  are 
extremely  useful  for  validating  numerical  fate  and  transport  mo^ls  which  include  complex 
processes,  such  as  solute  advection-dispersion,  geochemical  attenuation,  and  microbial 
degrulation.  A  secondary  objective  was  to  measure  the  in  situ  degradation  rates  of  the 
selected  organic  compouiids,  in  order  to  investigate  the  possibility  of  utilizing  natural 
attenuation  processes  as  an  alternative  to  active  remediation  techniques. 

BACKGROUND 

Groundwater  contamination  is  an  important  environmental  concern  for  both  military 
and  civilian  industrial  qierations.  The  organic  constituents  of  fuels  such  as  JP-4,  gasoline, 
and  diesel  (e.g.,  benzene,  toluene,  naphthalene,  xylene,  etc.)  and  cleaning  solvents  (e.g., 
trichlmo^^lene)  have  resulted  in  serious  corttamination  problems  at  many  DOD  facilities. 
Many  of  these  chemical  conqronents  of  fuels  and  solvents  are  regulated  by  EPA  as  hazardous 
substances,  and  have  stringent  concentration  limits  in  drinking  water  (Geraghty  and  Miller, 
Inc.,  1991).  Their  presence  in  groundwater  at  concentrations  exceeding  the  ^A*s  maximum 
contaminant  levels  generally  requires  the  use  of  costly  water  purification  methods.  The  lack 
of  sm  adequate  quamitative  understanding  of  the  processes  that  govern  the  environmental  fate 
of  these  contaminants  has  constrained  the  development  of  cost-effective  groundwater 
remediation  techniques. 

Degradation  nde  constants,  for  use  in  fue  and  transport  models  of  reactive  organic 
conqxHinds,  must  be  determined  by  experin^ntation.  Field  experiments  to  determine 
degnaiation  rates  are  desirable  because  laboratory  vali^  may  not  be  valid  for  the  actual 
condititHis  in  an  aquifer.  Madsen  (1991)  notes  diat  the  "determination  of  microbial  activity 
in  disturbed,  displaced  environmental  samples  incubated  in  the  labmatory  is  likely  to  be 
quantitatively,  even  qualitatively  different  from  the  same  determination  in  situ.  ” 

Determination  of  a  degradation  rate  in  an  aquifer  should  meet  criteria  that  are  iqiplied  to  rate 
measuranents  in  laboratory  batch  reactor  systems. 

Therefore,  a  known  mass  of  a  pure  conqxiund  must  be  introduced  into  the  aquifer  at  a 
^lecific  location  and  time,  followed  by  observation  of  the  changes  in  concentration  over 
^pace  and  time.  Because  aquifers  are  c^n  systems,  water  samples  should  be  analyzed  for 
both  the  reactive  compound  and  its  degnaiation  products,  and  a  mass  balance  should  be 
maintsuned  for  all  of  these  components.  Since  a  rate  constant  calculated  from  the 
disiqjpearance  of  the  contaminant  may  represent  physical  losses  rather  than  chemical 
dc^iradation,  the  demonstration  of  degradation  in  an  aquifer  requires  the  identification  and 
^landfication  of  reaction  products.  Field  measurements  of  organic  compound  degradation 
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rates  Att  meet  these  criteria  were  not  found  in  the  literature.  The  natural  attenuation 
eaqucriaacnt  conduced  by  Barker  et  al.  (1987)  measured  mass  loss  rates  of  aromatic 
oorapouods  from  a  pulse  plume,  but  did  not  establish  that  observed  l<»ses  were  due  to 
d^radatkm. 

An  eariier  natural  gradieitt  experiment,  referred  to  as  MADE-1,  involving  brcmiide 
and  three  fhKwobenzoate  tracers,  was  conducted  between  October  1986  and  June  1988  at  the 
same  site.  Boggs  et  al.  (1992)  described  the  test  site,  experimental  design  and  procedures, 
and  {wesent  cpialitative  observations  regarding  the  migration  behavior  of  the  bnnnide  tracer 
phone.  A  oonprehensive  description  of  the  hydrogeology  of  this  site  is  given  by  Boggs 
et  al.  (1990). 

Adams  and  Gelhar  (1992)  provided  a  spatial  moments  analysis  of  the  bromide  data  for 
die  MAl^-1  field  experiment.  Their  analysis  indicates  that  the  bromide  mass  recovery 
decreased  to  i^ifHXiximately  SO  percent  by  the  end  of  the  study  due  to  the  combined  effects  of 
iduroe  truncation,  sanqilii^  bias,  and  sorption.  However,  th^  showed  that  the  relationship 
qS  ho^itudinal  variance  and  mean  phune  displacement  was  relatively  insensitive  to  the 
bcmnide  mass  loss  during  the  experiment,  .^iplicatkm  of  a  nonuniform  advection-di^ierskin 
model  to  die  brcmiide  phune  moments  lead  to  an  estimate  of  longitudinal  di^rsivity  of  S  to 
10  mmers  for  die  alluvial  aquifer  (Adams  and  Gelhar,  1992). 

During  MADE-1,  ^  btnehole  flowmeter  was  demonstrated  to  be  a  practical  device 
fm  obtainmg  the  three-dimensional  distribution  of  horizontal  hydraulic  ccmductivity. 

Ririifiridt  et  al.  (1989a)  evaluated  the  accuracy  and  reliidiility  of  the  tedinkpie  through  a 
series  of  field  and  b^ratory  tests,  and  provkied  diemetical  and  practical  guidelines  for  its 
inoidementation.  .^pidkation  of  the  flowmeter  and  odier  methods  for  estimating  hydraulic 
conducdvi]^  variability  and  aquifer  diqiersivity  at  the  test  site  is  discussed  by  Rdtfeklt  et  al. 
(1989b.  1^). 

SCOPE 

The  locadon  of  this  field  study  at  the  Cohimbus  AFB,  Columbus,  MS, 
macrodi^NUsion  expCTiment  (MADE)  site,  is  shown  in  Figure  1.  In  this  natural  gradient 
experiment,  referred  to  as  MADE-2,  the  researchers  made  a  pulse  injection  of  a  dilute 
aqiieomt  whttkm  of  tritialed  water,  benzene,  p-xylene,  mqdtfhalene  and  o-dichlorobenzene 
(0-D^).  the  injection  was  made  into  the  satunbed  zone  of  the  unconfined  aquifer.  The 
dqpfndatkm  rates  of  diese  con^Kxmds  in  the  Columbus  aquifer  were  calculated,  and  these 
rates  were  related  to  the  aqaifo’s  structure  and  hydrologic  properties. 

The  MADE-2  experiment  was  initiated  w^  a  2-day  pulse  injection,  beginning  on 
Ame  26,  1990.  Tritiated  water  tracer  and  organic  solute  ccmcentration  disQibutions  were 
sidrseqiientty  mwiiored  in  three  dimensions  at  1-  to  3-month  intervals  over  15  mcmdis,  usii% 
im  exWBsive  netwmk  of  saturated  zone  multitevel  samplers.  Selected  sanqrlets  were 
monitmed  more  frequently  during  the  experiment  to  develop  complete  concentration  time- 
series  at  several  points  located  in  the  path  of  the  st^te  plume. 

Miration  bdiavior  of  the  solutes  during  tire  experiment  was  prinoarily  determined  by 
qratial  moosmits  Sfuitysis  of  tire  concetttration  measurements.  This  method  permitted 
coinpatnoii  of  tire  mi^  features  (e.g.,  relative  mass,  mean  di^lacement,  ^tial  covariance, 
dEebness,  md  kurtosis)  of  the  various  solute  plitines.  Estimates  of  the  relative  mass  (zero'" 
moment)  of  the  conservative  tracer  during  the  study  demonstrated  the  overall  effectiveness 
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Figure  1.  Test  Site  Location  Map. 
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of  die  sMo^ing  program.  Application  of  a  nonuniform  advection-dispersion  model  to  the 
fkst  diMe  laoinenls  of  the  conservative  tracer  data  was  used  to  estimate  the  apparem  aquifer 
tdapcwivhy  for  die  site.  CooqwristMis  between  the  relative  mass  estimated  for  die 
comervadve  inomr  tai  the  otpaic  sotules  fMOvided  a  means  (tf  quantifying  (xganic  mass 
losses  dne  to  the  condiined  effects  of  soiptkm  and  degradatkm. 

A  aepttiie  anafysis  of  dme-series  diaa  at  selected  points  was  done  to  estinuoe 
rettfdatioii  fectors  and  degradation  rates  for  the  organic  sohifes.  Field  estimates  of  sorpckm 
coefficients  for  diese  organic  ctunpounds.  derived  from  this  time-series  analysis,  were 
ctMBfMeed  widi  estiinates  obtained  from  labor»ory  batch  and  column  ejqierimeitts.  Results  of 
diese  congNtfisonB  he^  establish  the  effectiveness  of  labtmitoiy  smrptkm  measuranents  for 
predict  freld-scale  sorptitm  bdiavimr. 
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SECTION  II 
SITE  DESCRIPTION 

HydrogeolQgic  and  geostatistical  descriptions  of  the  aquifer  at  the  test  site  have  been 
presented  in  several  earlier  publications  (e.g.,  Boggs  et  al.,  1990  and  1992;  Rehfeldt  et  al., 
1992).  Relevant  results  from  these  works  are  summarized  below,  and  provide  the  reader 
background  information  for  use  in  interpreting  the  descriptions  of  solute  behavior  given  in 
this  report. 

AQUIFER  CHARACTERISTICS 

The  shallow  unconfined  aquifer  which  inunediately  underlies  the  site  consists  of 
alluvial  terrace  deposits  averaging  approximately  11  meters  in  thickness.  The  aquifer  is 
coiiq)osed  of  pooriy-  to  well-sorted  sandy  gravel  and  gravelly  sand  with  variable  silt  and  clay 
content.  Sediments  are  generally  unconsolidated,  and  occur  as  irregular  horizontal  or  nearly 
horizoidal  lenses  and  layers.  Marine  sediments  belonging  to  the  Eutaw  Formation  and 
consisting  of  clays,  silts,  and  fine-grained  sands,  form  an  aquitard  beneath  this  alluvial 
acpiifer.  Hydrogeologic  properties  of  the  Columbus  aquifer  and  a  desci^tion  of  the  MADE 
site  were  given  by  Boggs  et  al.  (1992). 

Aquifer  porosity,  bulk  density,  and  particle  density  were  determined  by  gravimetric 
and  vohimetiic  analyses  of  84  minimally  disturbed  soil  cores  by  Boggs  et  al.,  (1990).  The 
sanq>le  mean  bulk  density  of  the  84  sanq>les  was  1.77  g/cm’,  and  the  standard  deviation  of 
the  measurements  was  0.18  g/cm^  The  mean  and  standard  deviation  of  the  particle  density 
measurements  were  2.S7  and  0.07  g/cm^  reflectively.  The  mean  and  standard  deviation  of 
the  porosity  measurements  were  0.31  and  0.08,  reflectively.  A  porosity  of  0.3S  is  assumed 
for  the  fiatial  omments  analysis  presented  later  in  the  report,  bas^  on  an  estimation  of  cote 
compaction  used  in  analysis  of  the  MADE-1  data  (Adams  and  Gelhar,  1992).  Boggs  et  al. 
(1990)  also  presented  data  showii^  the  wide  range  of  aquifer  material  particle  sizes  and  the 
spatial  heterogeneity  of  the  particle  size  distrilwtions  within  the  Columbus  aquifer. 

Geostatistical  analysis  of  2187  measurements  of  the  horizontal  radial  component  of  the 
hydraulic  conductivity,  k],,  by  the  borehole  flowmeter  method  (Rehfeldt  et  al.,  1989a) 
indk»tes  that  die  variance  of  the  natural  logarithm  of  K,,  is  4.S  for  the  aquifer  (Boggs  et  al., 
1990;  Rehfeldt  et  al.,  1992).  Hereafter  in  this  report,  is  referred  to  as  the  hydraulic 
conducdviQr,  but  this  usage  refers  to  the  results  of  the  bore-hole  flow  meter  measurements, 
and  does  not  represent  a  ctmiplete  hydraulk;  conductivity  tensor.  Horizontal  and  vertical 
correlation  scales  of  12.8  and  1.6  meters,  ief)ectively,  were  also  estimated  from  these  data. 

The  Kh  profile  and  the  vertically-averaged  map  presented  on  Figure  2  illustrates 
the  extrrane  hetero^neity  of  this  aquifer.  Hydraulic  conductivities  measured  at  discrete 
vertKal  intervals  in  fully  screened  test  wells  typically  range  over  two  to  four  orders  of 
magnitude  at  each  test  well.  The  spatial  distri^tion  of  hydraulic  conductivity  varied  from 
sqqproximately  10^  to  K/*  cm/s.  The  magnitude  of  this  variability  is  apparent  from 
comparison  with  the  much  less  variable  (10^  to  10'^  cm/s)  field  at  the  Borden  Canadian 
Forces  Base,  Ontario,  described  by  Robin  et  al.  (1991).  The  mean  hydraulic  conductivity 
along  the  tracer  travel  path  increases  from  approximately  lO*’  cm/s  near  the  injection  point, 
or  near  field,  to  lO*^  cm/s  and  larger  in  the  far-field.  At  the  northern  extreme  of  the  test 
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Figttie  2.  (a)  Plan  View  of  V(ntiadly  Averaged  Hydraulic  Conductivity,  and 
Hydraulic  Conductivity  Di^bution  Along  Section  A-A’. 
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site,  the  mean  again  decreases  to  about  10'^  cm/s.  The  central  region  of  relatively  high 
mean  K|,  roughly  corresponds  to  a  former  river  meander.  This  meander  trends  southwest  to 
northeast,  and  is  at  about  a  4S-degree  angle  to  the  general  direction  of  groundwater  flow 
(Figure  3).  However,  tracer  tests  indicate  that  groundwater  flow  does  not  follow  the  axis  of 
the  meander  channel  as  might  be  expected  if  the  channel  represented  a  regionally  continuous 
zone  of  relatively  high  mean  Kh  (Rehfeldt  et  al.,  1992).  Note  that  all  elevations  presented  in 
this  report  are  referenced  to  mean  sea  level  (MSL). 

The  large-scale  spatial  variations  in  mean  K,,  at  the  test  site  are  related  to  the  region 
of  horizontally  converging  groundwater  flow  shown  on  Figure  3.  The  zone  of  convergent 
flow  is  associated  with  relatively  high  hydraulic  conductivity  sediments  present 
downgradient  of  the  injection  point.  Based  on  the  observed  migration  of  the  bromide  plume 
during  MADE-1,  groundwater  velocity  increases  downgradient  from  this  zone  of 
convergence.  The  mean  groundwater  velocity  probably  decreases  near  the  northern  end  of 
the  sampling  networic  due  to  a  decrease  in  mean  K^. 

SOIL  MINERALOGY 

Sand  and  gravel  particles  in  the  aquifer  are  composed  of  quartz,  feldspar,  and  mica. 
Fine-grained  materials  (i.e.,  those  passing  the  0.074  mm  sieve)  consist  of  quartz,  potassium 
feld^Nur,  muscovite,  and  clay  minerals.  Clay  minerals  present  are  predominantly  kaolinites 
and  Ulites  with  minor  amounts  of  montmorillonites  and  vermiculites.  Soil  particles  are 
commonly  coated  with  iron  oxides.  Measurements  of  free  iron  oxides  range  from  2.0  to 
3.9  percent  by  weight  for  particles  less  than  2  nun  in  diameter,  and  from  O.S  to  0.8  percent 
for  particles  larger  than  2  mm  (Boggs  et  cU.,  1990).  MacIntyre  et  ai.  (1991)  report  low 
organic  carbon  contents  for  SO  samples  of  alluvium  with  measurements  ranging  from  0.02  to 
0.06  percent  by  weight. 

GROUNDWATER  CHEMISTRY 

Some  chemical  characteristics  of  the  groundwater  from  wells  in  the  alluvial  aquifer  at 
the  site  prior  to  this  test  are  summarized  in  Table  1.  The  groundwater  had  a  low  total 
dissolved  solids  content  averaging  43  mg/L.  The  ionic  composition  of  the  groundwater  was 
^Mninat^  by  sodium,  silica  and  chloride.  Total  acidity  averaged  71  mg/L  (as  CaCOj). 
Alkalinity  averages  10  mg/L,  mainly  in  the  form  of  bicarbonate.  Dissolved  carbon  dioxide 
was  quite  high  at  the  groundwater  pH,  which  averages  4.8.  The  free  CO2  content  averaged 
57  mg/L,  representing  about  3  percent  of  the  total  dissolved  gas  concentration.  The  mean 
Qi  of  543  mV  indicates  an  oxidizing  groundwater  environment  (Boggs  et  al.,  1990).  These 
observations  of  groundwater  chemistry  are  preliminary  and  only  intended  to  describe  the 
general  chemical  conditions  the  solutes  might  encounter.  The  number  of  wells  sampled  was 
not  sufficient  to  allow  estimation  of  the  statistical  significance  of  the  measurements,  and  the 
set  of  chemical  parameters  sampled  did  not  allow  estimation  of  closeness  to  redox 
equilibrium  in  the  groundwater.  A  more  complete  geochemical  survey  of  the  site  is  now 
beii^  carried  out  to  support  micobiological  experiments  planned  for  a  new  test  at  the  site. 
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TABLE  1.  SUMMARY  OF  GROUNDWATER  CHEMICAL  CHARACTERISTICS 


mean 

min 

max 

N 

pH 

4.8 

4.3 

5.2 

21 

Eh  (mV) 

543 

410 

672 

16 

conductivity  (mmhos/cm) 

45 

25 

74 

16 

temperature  (°C) 

17.9 

9.8 

20.9 

18 

total  dissolved  solids  (mg/L) 

43 

30 

70 

6 

acidiQr  (mg/L  CaCOj) 

71 

36 

93 

11 

alkalinity  (mg/L  CaCOj) 

10 

2 

35 

15 

Ca-Mg  hardness  (mg/L  CaCOj) 

7.8 

7.7 

8.0 

2 

total  organic  carbon  (mg/L) 

2.3 

“ 

— 

1 

sodium  (mg/L) 

5.3 

— 

“ 

1 

silica  (mg/L) 

5.0 

— 

“ 

1 

calcium  (mg/L) 

1.8 

1.7 

1.8 

2 

magnesium  (mg/L) 

0.84 

0.77 

0.92 

2 

potassium  (mg/L) 

0.6 

— 

“ 

1 

iron  (Mg/L) 

15 

-- 

— 

1 

copper  (mg/L) 

7.5 

5 

10 

2 

zinc  (mg/L) 

2.6 

-- 

" 

1 

manganese  (mg/L) 

1.5 

— 

1 

chloride  (mg/L) 

8.0 

6.0 

10.0 

2 

sulfate  (mg/L) 

2.0 

1.0 

3.0 

2 

nitrate-nitrite  (mg/L) 

1.5 

1.1 

1.8 

2 

bromide  (mg/L) 

0.05 

0.04 

0.06 

7 

sulfide  (mg/L) 

8 

7 

8 

2 

total  dissolved  gas  (%) 

110 

105 

114 

6 

dissolved  oxygen  (mg/L) 

4.2 

1.6 

6.7 

14 

free  CO2  (mg/L) 

56 

56 

57 

2 

CO2  (mg/L) 

73 

58 

88 

2 
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SECTION  m 
METHODS 

TRACERS  AND  REACTIVE  SOLUTES 

Solutes  selected  for  the  experiment  are  listed  in  Table  2,  along  with  the  initial 
concratiation  and  total  injected  mass  (activity)  for  each  solute.  Tritiated  water  served  as  the 
conservative  reference  tracer.  The  four  organic  compounds  chosen,  benzene,  naphthalene, 
p>xylaie,  and  o^lichlorobenzene  (o-DCB),  are  common  constituents  of  various  hiels  and 
soivoits.  These  conqxxinds  vary  considerably  in  their  mobility  in  aquifer  materials  and  in 
their  suscq[>tibility  to  biotransformation  in  laboratory  systems.  Tabak  et  al.  (1981)  reported 
Mgnifirant  degradation  with  rapid  adaptation  for  laboratory  studies  involving  benzene  and 
n^q|Athalene,  wheicas  o-dichlorobenzene  showed  slower  rates  of  degradation  with  gradual 
adequation  followed  by  a  de-adaptive  process.  A  small  fraction  of  the  injected  p-xylene  was 
ring-hdieled  with  to  permit  measurement  of  the  degradation  of  p-xylene  to  its  reaction 
products. 

TABLE  2.  INITIAL  SOLUTE  CONCENTRATIONS  AND  MASSES  (ACTTVITIES) 

Mean  Mass/ Activity 


Tracer 

Concentrathm 

Ii^ected 

tritium 

55,610  'gCyvdL 

0.5387  Ci 

^*C  (p-xylene) 

2770  pCi/mL 

0.0268  Ci 

benzene 

68.1  mg/L 

659.7  g 

p-xylene 

51.5  mg/L 

402.0  g 

naphthalene 

7.23  mg/L 

70.0  g 

o-dichlorobenzene 

32.8  mg/L 

317.7  g 

SOLUTE  INJECTION 

The  rnttbod  of  solution  injection  was  similar  to  that  used  in  the  earlier  MADE-1 
brmnide  experimem  (Boggs  et  al.,  1992).  The  solution  was  injected  through  five  wells 
spaced  l-metor  q)art  in  a  linear  array  (Figure  4).  Each  injection  well  was  screened  over  a 
O.d-meter  interval  between  the  elevations  of  57.5  and  58. 1  meters.  The  entire  injection  was 
over  a  pmiod  of  48.5  hours  beginning  June  26,  1990,  at  a  uniform  rate  of  3.3 
f  /mipi  The  total  vohune  of  solution  injected  was  approximately  9.7  cubic  metms.  The 
maxinuim  increase  in  hydraulic  head  near  the  injection  wells  was  0.45  meters. 

The  iqjgctioff  strfittion  was  prepared  and  stored  on-site.  Ambient  groundwater  from  a 
wdl  lor^fd  apiuoximately  75  meters  upgradient  from  tlte  injection  point  was  used  to  prq>are 
die  The  trUium  and  pure-phase  organic  compounds  were  mixed  and  stored  in  two 

thonosiAy  2.7  m^  vcUume  taidcs,  each  equijq)^  with  a  floating  lid  and  an  impeller 

mixer.  Sanqdes  (rf  the  injection  solution  were  collected  at  1-  to  2-hour  intervals  from  the 
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Figure  4,  Map  of  Well  Network  Showing  Injection  Wells,  MLS,  and  BarCad*  Samplers. 
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8Hms  leadiBg  co  the  inyecttoa  weils.  The  injected  mass  (activity)  for  each  solute  was 
catealated  from  the  concentration  measurements  for  the  samples  taken  duriitg  injection. 

SOLUTE  M(WrFORING 

Figure  4  also  shows  the  locations  of  the  328  multilevel  samplers  (MLS)  used  to 
monitor  the  solute  plumes  in  three  dimensions  during  the  experiment.  Locations  of  sampling 
poiitts  are  defined  in  a  three-dimensional  Cartesian  coordinate  system,  with  the  origin  at  the 
ceitter  of  the  iiyection.  The  x  axis  is  approximately  normal  to  the  plume  trajectory,  the  y 
axis  is  tq^roximately  parallel  to  the  plume  trajectory,  and  the  z  axis  represents  the  elevation 
of  sampling  points  i^ve  MSL.  Specifically,  the  positive  y  axis  points  in  a  direction  12 
degrees  West  of  true  North. 

Details  regarding  the  design  and  constrw^n  of  the  mtiltilevel  samplers  are  given  by 
Boggs  et  al.y  1988.  Each  MLS  was  equipped  with  20  to  30  sampling  points  spaced  0.38 
m^ers  apart  in  the  vertical  dimension.  Groundwater  samples  from  each  MLS  were 
simultaneously  collected  using  a  mobile  sampliqg  cart  equii^)ed  with  three,  lO-chanoel 
per»taltic  pumps.  Sanqding  protocols  tetputed  the  puiging  of  120  mL  (approximately  four 
sample  tute  volumes)  of  groundwater  prior  to  die  collection  of  a  25  mL  sample  in  a  glass 
vial  with  a  Teilon*-lined  sqptum  cap.  Before  sealing,  0. 1  mL  of  a  10  percent  sodium  azide 
sohitiott  (a  hacteryostat)  was  added  to  each  water  sample.  Finally,  they  were  refrigerated  at 
tq)|»pximatdy  and  stored  in  an  inverted  petition  (i.e.,  capp^  end  down)  until  analyzed. 
Five,  diree-dimensioiial  sampling  events  (snapshots)  of  the  solute  plumes  were  performed  at 
apimiximatdiy  IQQ-day  intervals  during  the  s^y.  General  information  related  to  these 
sanqiliQg  events  is  given  in  Table  3.  Note  that  the  fifth  snapshot,  conducted  at  440  days, 
only  defined  the  organic  solute  plumes,  which  occurred  in  tte  region  near  the  injection  wells 
and  did  not  cross  the  well  field  boundaries.  At  440  days,  the  tritium  and  ‘*C  plumes  had 
crossed  the  downgradient  boundary  of  the  well  field,  and  thus  were  not  completely  sanqtled 
during  the  fifth  snapdiot. 
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TABLE  3.  GROUNDWATER  SAMPLING  SUMMARY 


Sampling 

Event* 

Date 

Elapsed 

Time  (days) 

No.  Wells 
Sampled 

Samples 

Analyzed 

FOl 

Jul9-ll,  1990 

13 

26 

264 

S21 

Jul  23-27,  1990 

27 

99 

1226 

F02 

Aug  13-17,  1990 

48 

31 

287 

F03 

Sep  17-19,  1990 

83 

53 

185 

F04 

Oct  5-17,  1990 

111 

39 

72 

S22 

Nov  5-8,  1990 

132 

111 

871 

F05 

Dec  3-4,  1990 

160 

29 

195 

F06 

Jan  8-9,  1991 

195 

25 

196 

S23 

Feb  5-7,  1991 

224 

190 

1976 

FO? 

Apr  3-5,  1991 

281 

42 

387 

S24 

May  21-23,  1991 

328 

205 

2345 

S25 

Sep  9-11,  1991 

440 

79 

460 

F  =  fenceline  event;  S  =  snapshot 

Fifty-six  positive  displacement  BaiCad*  water  samplers,  manufactured  by  Goldberg, 
Zoino  and  Associates,  Newton  MA,  were  installed  in  wells  that  formed  two  parallel  rows  or 
"fencelines”  oriented  normal  to  the  genera!  flow  direction.  Their  locations  are  also  shown  in 
Figure  4.  The  BaiCad*  sanq>lers  were  permanently  installed  at  fixed  depths  in  uncased, 
back  filled  boreholes.  Several  samplers  were  installed  in  each  borehole  to  allow  sampling  at 
multiple  depths. 

Tlw  BarCad*  sampler  contains  a  check  valve  and  porous  filter  through  which  water 
passes  before  exiting  to  die  surface  in  a  small  riser  tube,  which  is  inside  of  a  larger  gas  drive 
tube  (0.5  and  1.3  cm  diameter,  respectively).  The  device  uses  a  closed-system  pressure 
drive  principle  (positive  inert  gas  displacement)  to  collect  the  sample  and  maintain  sample 
integrity.  Ceramic  filters  were  selected  for  this  aj^lication  based  on  the  results  of 
laboratory  tests  of  sampler  efficiency,  using  standi  hydrocarbon  solutions. 

The  two  fence  rows  were  located  at  distances  of  approximately  6  and  16  meters  from 
the  injection  point.  The  vertical  distribution  of  the  sampling  points  across  each  row  is  shown 
in  Figure  5.  BarCad*  sampling  devices  were  installed  in  IS  cm  diameter  auger  holes  with  no 
more  than  two  samplers  per  hole.  Bentonite  seals  approximately  0.6  meters  thick  were 
placed  in  the  annuls  space  above  and  below  each  sampling  point  to  hydraulically  isolate 
sampling  points  install^  in  the  same  auger  hole. 

Sandies  from  the  fence  rows  of  BarCad*  points  were  taken  at  frequent  intervals  to 
provide  time-series  data  for  analysis  of  organic  solute  sorption,  retardation  and  degradation. 
The  sampling  frequency  varied  from  2  week  intervals  initially  to  3  month  intervals  later  in 
the  study,  as  shown  in  Table  3.  Time-series  samples  were  also  collected  at  a  few  MLS 
poirtts  concurrently  with  the  BarCad*  fence  row  samplings.  Results  from  both  of  the 
sanqilers  were  used  to  determine  the  similarity  of  hydrocarbon  sorption  and  degradation 
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panmeters  calculated  from  time-series  MLS  data  and  from  time-series  BaiCad*  data.  This 
test  of  sampler  similarity  is  described  in  more  detail  in  Section  6. 

As  shown  in  Figure  4,  the  spatial  density  of  MLS  wells  near  the  injection  wells  is 
high,  and  BarCad*  sample  points  in  the  fence  row  wells  are  close  to  MLS  sample  points  in 
die  neariiy  MLS  wells.  This  proximity  enables  comparisons  to  be  made  between 
concentrations  of  solutes  obtained  simultaneously  by  the  MLS  and  BaiCad*  samplers  during 
the  field  experiment.  The  BaiCad*  samplers  thus  provided  an  indirect  test  of  die  probability 
of  solute  loss  during  sampling  due  to  sorpdon  and  vaporization.  Since  the  principles  of 
operation  and  materials  of  construction  differed  between  the  two  sampler  types,  the  observed 
agreement  of  concentration  results  for  samples  from  proximal  BaiCad*  and  MLS  sample 
points  indicates  that  losses  associated  with  a  particular  sampler  were  probably  small.  If  such 
losses  did  occur,  it  is  unlikely  that  they  would  be  equal  for  the  different  sample  devices.  The 
use  of  BarCad*  saiiq)lers  ratter  than  MLS  samplers  in  the  fence  rows  served  three  purposes: 
(1)  to  demonstrate  the  utili^  of  BaiCad*  samplers,  (2)  to  obtain  time-series  concentration 
data,  and  (3)  to  demonstrate  that  solute  losses  during  sampling  did  not  pose  a  problem.  The 
third  purpose  is  obviously  critical  to  attainment  of  the  first  two,  and  fortunately,  the  BaiCad* 
data  compared  well  with  MLS  data  from  adjacent  wells.  This  agreement  could  be  confirmed 
by  a  nearest  neighbor  con^iarison  of  solute  concentrations  between  BaiCad*  and  MLS  well 
points  ,  but  this  task  would  be  extremely  tedimis.  An  alternative  comparison  is  made  here, 
based  on  die  presence  (see  Fipire  4)  of  rows  of  MLS  wells  immediately  upgradient  and 
down  gradient  of  the  BarCad*^ fence  rows. 

At  each  sampling  time,  an  average  coiteentration  for  each  solute  was  computed 
separately  for  all  of  the  BaiCad*  sampling  points  in  each  fence  row,  and  for  the  sampling 
points  in  the  associated  upgradient  and  downgradient  MLS  rows.  Results  for  the  fence  row 
at  6  meters  from  the  injection  wells  are  given  in  Table  4.  Comparison  of  the  average 
concentrations  of  solutes  in  the  two  sanqiler  types  shows  that  the  results  are  reasonable,  with 
a  general  trend  of  decreasing  values  fn»n  the  upgradient  MLS  row  to  the  BaiCad*  fence  row, 
and  finally  to  the  downgradient  MLS  row.  Th^  findings  afqiear  to  support  the  utility  of 
BarCad*  sanqilers  for  such  a  purpose.  Unforhinately,  this  comparison  must  be  mainly 
subjective  due  to  the  heterogeneity  of  the  aquifer  and  the  unavoidable  spatial  separation  of 
MLS  and  BaiCad*  sampling  points. 

ANALYTICAL  METHODS 

Stable  organic  compounds  in  water  samples  were  analyzed  by  gas  chromatography 
(GC).  D^terated  toluene  and  4-bromofluorobenzene  internal  standards  were  added  to  20 
mL  water  samples,  which  were  then  extracted  with  2  mL  of  n-pentane.  The  extracts  were 
analyzed  for  benzene,  p-xylene,  naphthalene  and  o-dichlorobenzene  by  (}C,  using  a  30-meter 
long,  0.32  mm  i.d.  fused  silica  column  with  1  ^m  thick  DB-S  bonded  phase  coating  and 
flame  ionization  detection.  Helium  carrier  gas  flow  was  set  at  3-S  mL/min  and  the  split  ratio 
was  10:1.  The  initial  column  temperature  was  held  at  50°C  for  1  minute  and  was  then 
increased  at  10*’C/min  to  ISO^C,  which  was  held  for  S.S  minutes.  Sensitivity  of  the  method 
was  50  ng/L  for  benzene  and  4  ng/L  for  the  other  compounds.  All  samples  were  extracted 
and  analyzed  within  21  days  of  collection. 

Quality  control  for  the  organic  analyses  was  ensured  by  using  several  techniques. 

First,  analytic  precision  was  checked  through  duplicate  analysis  of  sample  extract  splits  for 
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4.  AVERAGE  SOLUTE  CONCENTRATIONS  IN  6  METER  FENCE  ROW 
BARCAD*  SAMPLERS  AND  IN  ADJACENT  ROWS  OF  MLS  SAMPLERS 
Oif/L  for  Hydrocarbons  and  pCi/mL  for  Radioisotopes) 
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[Napl 

[p-Xyl] 

[o-DCB] 

pH] 

V*C] 

83 

416.0 

26.2 

226.2 

176.8 

611.4 

6.1 

232 

573.6 

57.4 

458.9 

406.3 

915.2 

34.0 

}60 

525.4 

73.6 

453.4 

468.7 

957.4 

35.5 

224 

387.5 

55.7 

271.5 

353.4 

701.1 

25.5 

279 
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41.2 

59.0 

268.5 

526.8 

20.5 

328 
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16  3 
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344.5 

13.2 
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14.8 
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9.1 
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pH] 

P^ 

83 

388.4 

24.1 
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16.9 

|32 

167.2 

16.1 

127.4 

153.8 

457.9 

12.3 

160 

122.2 

14.0 

81.9 

133.7 

349.2 

9.9 

224 

191.9 

18.1 

123.6 

154.7 

400.1 

11.5 

279 

121.2 

12.3 

60.2 

113.1 

295.6 

9.3 

328 

61.7 

4.6 

7.5 

57.4 

486.0 

15.5 

440 

86.5 

11.6 

18.3 

90.2 

254.0 

7.3 

bowninadimn  MLS 

row 

IXdays) 

[Ben] 

[Nap] 

[p-Xyl] 

[o-DCB] 

pH] 

P^ 

tn 

235.4 

25.9 

163.3 

197.2 

383.4 

14.7 

224 

127.3 

14.5 

62.3 

97.0 

230.1 

9.8 

328 

81.3 

8.1 

11.9 

46.0 

158.0 

6.9 

440 

104.8 

10.7 

13.4 

77.4 

193.4 

8.2 

S  perorat  of  die  field  sanqiles.  Second,  certified  standard  solutions  obtained  fnmi  an 
iadt^endent  source  were  used  to  make  100  ftg/L  control  standards  for  each  organic  solute. 
These  eontred  ttandards  were  analyzed  at  a  frequency  erpiivalent  to  approximately  S  percent 
of  die  field  sanqiles  to  verify  analytical  accuracy.  Third,  d-bromoflourobenzene  was  added 
Ml  Ml  internal  standard  fen'  all  analyses,  and  deuterium-labeled  toluene  was  used  as  a 
Sfnrto|Bie.  Finally,  dedicate  field  samples  represetting  approximately  5  percent  of  the  total 
ndd  sanqdes  collected  for  each  sampling  evett  were  analyzed  by  the  Mississippi  State 
Chendcsl  UdKntttny  located  at  Mississ^i  State  University  (MSU),  as  an  indqiendent 
verifieitimi  of  the  Mudytkal  results. 

Trkhim  and  in  samples  from  the  wells  were  analyzed  by  liquid  scintillation 
ommting  in  dual-isoiope  mode,  by  die  Water  Resources  Research  Center  at  MSU.  A  5.0  mL 
siiqiiiot  of  each  sample  was  emulsified  in  15  mL  of  Pico-Flour  40  (Packard  In^rument  Co.) 
oodtttll.  Mid  hdd  in  darkness  for  at  least  1  hour  prior  to  counting.  Sanqiles  were  counted  on 
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a  Packard  Model  2250A  for  20  minutes  or  to  a  1  percent  error  at  the  95  percent  probability 
level,  whkhever  was  attained  first.  Background  levels  of  tritium  and  in  ambient 
groundwater  at  the  test  site  dictated  the  analytical  sensitivity  for  these  experimental 
measurements,  and  were  approximately  2  and  3  pCi/mL,  respectively.  The  statistical  errors 
associated  with  20-minute  measurements  at  these  levels  are  about  ±20  percent  at  the  95 
perceitt  [nobability  level. 

A  set  of  tritium  and  standards,  including  one  National  Bureau  of  Standards 
traceable  standard,  was  counted  every  108  samples  to  check  aiudytical  ruxuracy  and  verify 
instrument  performance.  Every  twenty-fifth  sample  was  counted  twice  to  estimate 
measurement  precision  and  repeatability.  As  an  overall  check  on  the  radiological 
measurements,  duplicate  field  samples  amounting  to  roughly  5  percent  of  the  total  field 
sanqil^  were  analyzed  independently  by  the  TV  A  Western  Area  Radiological  Laboratory. 

The  quality  control  data  for  the  radiological  and  stable  compound  analyses  can  be  obtaiiied 
from  the  Tennessee  Valley  Authority  Engineering  Laboratory,  Norris  TN. 

DISSOLVED  OXYGEN  MONITORING 

Dissolved  oxygen  (DO)  measurements  were  performed  at  2-  to  3-month  intervals  at 
selected  locatioiis  during  the  eiqKriniem.  DO  data  were  intended  to  demonstrate  whether 
aerobic  conditioiis  for  possible  degradation  of  the  organic  conqiounds  were  maintained.  As 
shown  in  Figure  6,  DO  monitoring  wells  included  16  MLS  located  along  the  longitudinal  axis 
of  the  tracer  plume,  and  6  piezmneters.  Three  DO  sanq>les  were  collected  from  each  MLS 
using  a  peristaltic  puiiq>.  The  mean  elevations  of  these  sample  points  were  iq>proximately 
S4.0,  57.4,  and  61.2  meters  above  MSL,  corresponding  to  the  lower,  middle  and  upper 
sections  of  the  aquifer.  Samples  were  collected  in  60  mL  BOD  bottles  after  purging 
iq)proximately  1()0  mL  from  the  MLS  tubes.  Measurements  were  made  with  a  calibrated  DO 
pr^  immediately  after  sample  collection.  The  DO  measurements  at  the  6  piezometers  were 
made  in  conjunction  with  pH,  conductivity,  and  oxidation-reduction  potential  measurements. 
Samples  were  collected  from  the  piezometers  using  a  bladder  pump,  and  DO  was  measured 
with  a  calibrated  in-line  probe.  Field  DO  data  and  coordinate  locations  of  the  DO  sampling 
points  are  given  in  Appendix  A. 
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SECTION  IV 

EXPERIMENTAL  OBSERVATIONS 


The  observed  transit  and  fate  of  tritium  and  organic  solutes  during  the  experiment 
are  qualitatively  represemed  and  ctescribed  in  this  section.  Evolution  of  the  various  solute 
phunes  in  q)ace  and  time  are  shown  in  a  series  of  graphs  of  the  spatial  concentration 
distributions  of  each  solute  for  each  snapshtH.  To  facilitate  direct  comparison  of  the 
transmit  behavior  of  the  various  solutes,  all  concentration  data  are  presented  in  terms  of 
relative  concentration,  i.e.,  concentration  have  been  normalized  to  their  respective  initial 
injection  concentration. 

Relative  (normalized)  concentration  are  generally  used  in  this  report.  The  raw 
concentration  and  other  field  data  can  be  obtained  by  request  from  the  Teiuiessee  Valley 
Authority  Engineering  Laboratory,  Norris,  TN. 

TRmUM  PLUME 

Vertical  profiles  (section)  of  tritium  activity  along  the  longitudinal  axis  of  the  plume 
for  the  first  four  siuqishots  are  presented  in  Figure  7.  Note  for  reference  that  ground  surface 
level  is  at  approximately  65  meters  above  MSL  in  diese  profiles.  A  similar  set  of  plots 
^wing  horizontal  section  through  the  tritium  phime  at  elevation  S9.S  meters  is  given  in 
Figure  8.  The  dmninant  feature  indicated  by  both  figures  is  the  extreme  asymmetry  of  the 
tritium  activity  distribution  in  the  longitudii^  dimenion.  In  this  respect,  the  tritium  plume 
was  quite  similar  to  the  MADE-1  bromide  plume  described  by  Adams  and  Gelhar  (1992). 
After  328  days,  the  most  concentrated  region  of  the  plume  remained  within  roughly  20 
meters  of  the  iiyection  point,  while  the  leading  edge  of  the  plume  extended  downgradient  a 
distance  in  excess  of  260  meters. 

This  tritium  migration  behavim  is  conistent  with  the  hydraulic  conductivity 
distribution  akmg  dte  friume  travel  path,  shown  in  Figure  9.  Tlie  slow  rate  of  displacement 
of  the  main  body  of  the  plume  in  the  near-Oetd  corresponds  to  the  relatively  low  mean 
hydraulic  conductivity  (aj^roximately  lO*’  cm/s)  in  this  region.  At  about  20-30  meters 
downgradieitt  of  the  injection  point,  the  mean  conductivity  of  the  aquifer  increases  by  one  to 
two  orders  of  magniUKte.  Based  on  the  260  meter  minimum  displacement  of  the  tritium  from 
within  224  days  of  injection,  groundwater  velocities  within  the  more  permeable  sediments  in 
the  far-fleld  in  excess  of  400  m/yr  are  indicated.  The  most  likely  conclusion  is  that  t}» 
extreme  skewness  of  the  tritium  plume  was  produced  by  accelerating  groundwater  flow, 
which  was  itself  caused  by  large-scale  qimial  variations  in  mean  hydraulic  conductivity. 

ORGANIC  SOLUTE  PLUMES 

Longitudinal  profiles  through  the  benzene,  naphthalene,  o-dichlorobenzene,  p-xylene, 
and  ^*C  plumes  for  each  snapshot  are  presented  in  Figures  10-14.  Corresponding  tritium 
l^ume  {uofiles  for  mich  date  are  given  at  the  top  left  of  each  figure  for  comparison  of  solute 
tranqiort.  Benzene  and  nafriithalene  showed  similar  transport  behavior  during  tlte  study. 

Both  had  limited  horizontal  displacemem  compared  with  tritium,  and  the  maximum  down- 
gradimit  extent  of  both  plumes  was  only  about  25  meters.  The  gradual  disappearance  of 
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Figiire  7.  Trithun  Relative<Activity  Profiles  in  the  Longitudinal  Direction. 
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Horizontal  Sections  Through  Tritium  Plume  at  Elevation  59.5  m  at  27,  132,  224,  and  328  Days. 
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Figure  9.  CQiiq>arison  of  Tritium  aixl  Hydraulic  ConductiviQr  Profiles 
in  the  Longitudinal  Direction. 
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Longitudinal  Profiles  of  Relative  Concentration  for  All  Solutes  at  132  Days. 
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Figure  12.  Longitudinal  Profiles  of  Relative  Concentration  for  All  Solutes  at  224  Days. 
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Figure  14.  Longitudinal  Profiles  of  Relative  Concentration  for  All  Solutes  at  440  Days. 


benzene  and  naphthalene  suggests  that  they  were  either  transformed  or  removed  from  solution 
dHring  Oe  experiment.  The  transport  behavior  of  o-DCB  was  similar  to  that  of  tritium  for  the 
first  132  days  of  die  field  study,  ^jor  differences  between  the  o-E>CB  and  tritium  plumes  were 
ftnt  evidcpt  at  224  days.  At  this  point,  the  o-DCB  plume  showed  less  overall  displacenMnt  and 
9atial  ctmdnuty  dian  the  tritium  plume. 

The  use  of  a  radiolabeled  compound  was  critical  towards  establishing  the  degradation 
of  these  test  chemicals.  During  the  study,  distinct  difierences  were  evident  between  the  p-xylene 
and  ^  plumes,  even  thmigh  the  **C  was  chemically  associated  with  the  injected  p-xylene.  The 
plume  diape  and  size  were  similar  to  that  of  the  tritium  plume  througlKMit  the  study.  In 
comiast,  p-xylmie  exhibited  limited  displacement  (less  than  SO  meters  in  the  down-gradient 
direction),  aid  diminished  plume  volume  and  mass  after  132  days.  Again,  transformation  of  the 
p-xylene  can  be  inferred  from  these  results. 

DISSOLVED  OXYGEN 

Field  dissolved  oxygen  (DO)  data  and  coordinate  locations  of  the  points  sampled  during 
the  experiment  are  given  in  Ai^ndix  A.  Temporal  average  DO  concentrations  in  the  region 
near  tte  iiyection  wells  for  each  well  sample  point  were  above  3.8  mg/L,  and  the  minimum 
uflividual  value  was  2.6  mg/L.  There  was  no  significant  difference  between  dissolved  oxygen 
cbnoratrations  bodi  near  and  far  from  the  injection  wells,  or  between  measurements  taken  before 
and  efier  infection.  This  result  was  antic^ated,  because  aromatic  hydrocarbon  concentrations 
in  die  aquifer  after  die  injection  were  too  low  to  significantly  dqilete  dissolved  oxygen  in  the 
aquifer. 

Figure  IS  ^ws  the  vertically  averaged  dissolved  oxygen  data  for  wells  along  the 
toagitudiiial  axis  of  the  tracer  plume  at  six  times  during  the  experhnent.  The  data  indicate 
conidendlle  tenqxnai  and  spatial  variability  in  DO.  Mean  values  ranged  from  4.3  to  6.0  mg/L 
during  die  study  period,  and  varied  ^latially  by  tqi  to  3  mg/L  over  this  transect  for  individual 
at^^ing  dates.  The  DO  near  the  iiqection  site  is  consisteiuly  lower  than  the  mean  DO  for  the 
transect.  Tliis  trend  is  evident  in  measurements  made  both  before  and  after  die  injection,  and 
nmy  nika  die  reladvely  slow  rate  of  groundwater  movement  in  this  region. 

CMssolved  oxygen  measurements  made  during  die  field  study  did  not  reveal  conclusive 
evideare  of  the  biological  vqnaiw  of  oxygen  within  the  organic  solute  plumes.  DO  time-series 
4$U  fitxn  MLS  points  within  the  more  concentrated  region  of  the  plumes  appear  to  indicate 
slight  DO  rechictions  of  tdiout  1  mg/L.  These  reductions  correlate  with  the  breakthrough  of  the 
mgaaic  solutes,  but  no  relation^p  can  be  proven  due  to  excessive  noise  in  the  DO  data, 
itc^fkms  of  greatm*  DO  uptake  may  have  exis^  near  the  injection  wells  in  the  first  few  day^, 
but  tliete  are  no  early  dam  to  test  this  possibility.  The  lack  of  significant  DO  uptake  may  also 
he  e]i|ihuned  by  eidianced  mixing  of  DO-rich  groundwater,  produced  by  the  extreme 
heterc^^eneity  of  die  aquifer.  The  importance  of  DO  transfer  by  dispersive  mixing  during 
iKKitranrfomiation  of  (»ganic  contaminants  has  been  demonstrated  in  the  numerical  experiments 
of  MacQuarrie  and  Sudicl^  (1990). 

There  is  no  evidence  that  anaerobic  conditions  occurred  at  any  time  in  this  experiment. 
Hie  geochemical  data  presented  earlier  indicate  that  the  nitrate  and  sulfate  concentrations  in  the 
mhb^nt  groundwater  are  insufficient  to  support  anaerobic  degradation  of  aromatic  hydrocarbons 
m  die  idume.  Groundwater  Eh  and  DO  met^urements  indicate  that,  throughout  the  experiment, 
the  a^ifer  environment  was  too  oxidizing  for  anaerobic  degradation. 
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Figure  IS.  Vertically  Averaged  Dissolved  Oxygen  Concentrations  at  Monitoring 
Points  Located  Along  the  Longitudinal  Axis  of  Solute  Plumes. 


As  a  check  m  the  reasonableness  of  the  fMld  DO  observations,  the  biological  oxygen 
dipmand  (BM>)  associirted  with  the  observed  aromatic  hydrocarbon  mass  losses  during  the  study 
was  csthnated.  The  reaction  stoichiometry  given  in  Equations  (1)  to  (4)  was  a5Maiin^  for  these 
odenh^tMis. 


benzene: 

2C^  +  1502  -  I2CO2  +  6H2O 

(1) 

p-xylene: 

2CjH,o  +  2IO2  -  I6CO2  +  IOH2O 

(2) 

npMialene; 

CioHg  +  I2O2  -  IOCO2  +  4H2O 

(3) 

o-DCB: 

CjH«Cl2  +  6.5O2  -  6CO2  +  2H^  +  2Cr  +  H2O 

(4) 

Estimates  of  total  BOD  for  each  snapshot  are  presented  in  Table  5.  To  c^tain  BOD  data 
in  concentntkm  units  for  comparison  with  field  DO  measurements,  it  was  necessary  to  Hrst 
estimate  the  volume  of  oxygenated  groundwater  that  had  mixed  with  the  organic  tracers  at  each 
sanqiling  dale,  then  use  Ais  volume  to  convert  oxygen  mass  to  oxygen  concentration.  For 
example,  initud  mixing  at  27  days  was  calculated  by  dividii^  the  tritium  plume  volume  of 
qi|»oxBnately  6  x  10^  liters,  by  the  injection  solution  volume  of  9700  liters,  to  obtain  a  dilution 
(rf  about  60. 


TABLE  5.  ESTIMATED)  BOD  REQUIRED  TO  M^RADE  ORGANIC  COMPOUNDS 


Elapsed  Tirae 

BOD 

Tl  nume 

(days) 

(mg/L) 

Estfanated  hfixing  Vol.  (L) 

27 

0.92 

0.5923  E+6 

132 

1.35 

0.1328  E-l-7 

224 

0.17 

0.1581  E+8 

328 

0.14 

0.2410  E-1-8 

The  volume  of  groundwater  associated  with  the  tritium  plume  for  each  snapshot  date 
was  assumed  to  be  approximately  equal  to  the  total  volume  of  groundwater  coming  in  contact 
witii  the  organic  solutes  for  purposes  of  the  calculation.  Results  indicate  that  His  volume- 
averaged  BOD  required  to  transform  all  of  the  oiganic  compounds  at  the  observed 
diisqipearance  rates  is  on  the  order  of  1  mg/L  or  less.  These  estimates  of  low  BOD  may 
^lain  the  nondetectability  of  DO  uptake  from  the  field  measurements. 
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HYDROLOGICAL  MEASUREMENTS 

The  maximum  seasonal  fluctuation  in  the  water  table  during  the  study  period  was 
approximately  3  meters,  as  shown  in  Figure  16.  This  change  resulted  in  about  a  30  percent 
variation  in  the  saturated  thickness  of  the  aquifer.  Mathematical  models  of  solute  transport 
must  therefore  include  a  time-varying  water  table,  since  it  is  generally  accepted  that  the 
variation  should  be  less  than  10  percent  for  application  of  the  assumption  of  a  time-invariant 
water  table.  Figure  16  also  presents  temporal  variations  in  the  magnitude  of  the  horizontal 
hydraulic  gradient,  |  J^l ,  the  direction  of  the  horizontal  gradient,  O^,  and  the  magnitude  of 
the  vertical  hydraulic  ^adient,  |  .  {1^1  was  estimated  by  fitting  a  plane  to  the  piezometric 

head  measurements,  and  averaged  0.0017,  with  seasonal  periodicity  corresponding  to  water 
table  fluctuations.  The  mean  and  standard  deviation  of  6,,  were  -4  degrees 
(cminterclodcwise  flom  the  y-axis)  and  28  degrees,  respectively.  The  general  direction  of 
the  hydraulic  gradient  is  consistent  with  the  observed  direction  of  the  tritium  plume 
migration,  as  shown  in  Figure  8.  The  relatively  large  transverse  dispersion  of  the  tritium 
plume  indicated  on  this  figure  may  have  been  partially  due  to  temporal  variations  in  the 
horizontal  direction  of  the  hydraulic  gradion.  The  vertical  hydraulic  gradients  are  equal  to 
or  greater  than  the  horizontal  gradients  at  several  locations,  and  are  generally  correlated  with 
water  table  fluctuations. 

The  temporal  changes  in  the  vertical  hydraulic  gradients  estimated  for  three  nested 
piezometers,  were  found  to  oe  dependent  on  their  location  at  the  field  site.  Typically, 
piezometers  located  upgradient  of  the  injection  point  exhibited  downward  (positive)  vertical 
grtulients.  The  magnitude  of  the  gradient  in  these  piezometers  occasionally  exceeded  10 
percem.  Vertical  gradients  indicated  for  piezometers  positioned  downgradient  of  the 
iiyection  point  generally  exhibited  a  lower  magnitude  and  mixed  directional  components. 

Tte  meaaired  vertical  gradients  during  the  experiment  were  generally  consistent  with 
observations  of  the  upward  movement  of  the  tritium  plume  near  the  injection  wells,  extensive 
vertical  mixing  farther  downgradient,  and  upward  motion  of  the  plume  near  the  downgradient 
boundary  of  the  well  field.  These  vertical  gradients  are  thus  in  agreement  with  the  tritium 
mass  transport  shown  in  Figure  7. 
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(a) 


(b) 


Figuie  16.  Temporal  Variatioiis  of  (a)  Potentioinetric  Surface  Elevation, 

(b)  M^^ntoKle  of  the  Horizontal  Hydraulic  Gradient  (  |Jh|  ), 

(c)  Directkm  of  the  Horizontal  Gndient  (  G|, ),  and 

(d)  Magnitude  of  the  Vertk»l  Hydraulic  Gradient  (  |Ji|  ). 
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SECTION  V 

SPATIAL  M(NMENTS  OF  SOLUTE  PLUMES 


The  spatial  moments  of  the  solute  concentration  measurements  provide  a  means  of 
characterizii^  and  ccmparing  the  tritium  and  organic  solute  plumes.  Inqwrtant  parameters 
determined  by  i^nttial  moments  analysis  include  the  total  solute  mass  (total  activity  of  tritium 
and  in  solution,  the  mean  plume  velocity,  and  the  spatial  covariance.  Because  the 
(dumes  had  strongly  non-Gaussian  shiq)es,  moments  calculations  of  the  skewness  and  kurtosis 
of  die  tracer  conceiUration  distributions  at  each  sample  time  were  done  to  further  defme  the 
shi^  of  CKb  phune. 

SPATIAL  MOMENTS  ESTIMATION 

The  formulation  of  the  spatial  mmnents  estimators  used  in  this  study  is  essentially  the 
same  as  that  enqiloyed  by  Adams  and  Gelhar  (1992)  in  their  analysis  of  the  MADE-1 
brtnnide  tracer  experiment.  Similar  methods  were  also  ad<^>ted  in  interpreting  other  recent 
field  tracer  eiqieriments  conducted  at  Borden  (Freyberg,  1986)  and  Cod  (Garabedian  et 
al.,  1991).  In  all  cases,  the  method  develcqied  by  Aris  (1956)  was  us^  to  relate  the  spatial 
moments  of  the  tracer  concentration  data  to  the  l^e-scale  plume  features.  The  general 
three-dimensional  spatial  moments  are  defined  by: 

“If!  ncxyzMxdydz,  (  x,  y,  z  -  oo  to  +  oe  )  (5) 

where  c  «  c(x,y,z,t)  is  the  tracer  concentration  roeasuremeitt,  n  is  porosity  (assumed 
constant  and  eq^  to  0.35  in  our  analysis),  and  i,  j,  and  k  are  imegers  ranging  from  zero  to 
infinily.  The  total  mass  associated  with  the  tracor  plume  is  given  by  the  zero'*  moment  about 


the  origin  (x  =  0,  y  0,  z  =  0)  is  given  by: 

Mono  =  1  1  J  ncdxdydz,  (  x,  y,  z  -  oo  to  +  oo  )  (6) 

The  first  inmnent  about  the  origin  provides  the  plume’s  center  of  mass  coordinates, 

Xc=  Mioq/Mqoo;  ye=  Moio/Mooo;  Moqi/Mqoo  (7) 

The  mean  plume  veloci^  vector  is  computed  fnnn  the  rate  of  translation  of  the  plume’s 
center  of  mass  (centroid)  as  shown  below, 

W,=  dxjdv,  V,=  dy,/dt;  W  =  dzjdt  (8) 

and  the  magnitude  of  the  mean  velocity  is  given  by: 

V  =  (V,HVy2+V/)®^  (9) 
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Ilie  seccnd  oicimeiit  about  the  coiier  of  mass  coordinates  of  die  plume  yields  the  following 
iqdBBfietric  spatial  covariance  tensor: 


'»  = 


®jof  ^xy 
Oyx  Oyy 


wbeit. 


M—  Mjim/Mooo  - 

o'^-NWMo.o-y.^ 

M002/M000- 


(10) 


xy“  ^ yx~~  Muo/Mfloo  *  X^c 
jj=  M|oi/Mooo*  x^ 

«^y*=  0*1  y*  Mo„/Mooo- 

Thfe  soMMid  immieitts  <^21.  ^  0^33)  ^  defined  1^  the  eigenvalues  of 

filiialkm  (10).  Tlie  skewness  and  Icurtmis  estiinatras  given  in  C^RC  (19S9,  p.  391)  were 
aSopK^  jibe  prmoit  anatysis.  Only  die  longitudinal  skewness  and  kurtosis  were 
eatimMtod.  iW  ineasine  of  longihkBnal  skewness,  g„  is  given  by: 


b\  “  (M(B(/Mooo  -  3MloioMoa(/M(xjo^  +  2Moio*/Mooo^)/2a’yy  (11) 

wfeiem  the  aimieild(»r  in  Equation  (11)  rqnexms  die  kmgitudinal  third  central  mcmient,  and  y 
crackles  wil&  die  longitudinal  axis  of  die  phime.  Longitudinal  kurtosis,  gj,  is  defined  by: 


*  (Mo«(M«»'^1^i<M>3o(®^ooo* +dMojo*Mo2o(Mooo*"3Mojol* /Mooo*)/2o^ yy"l  (12) 

wlkrib  the  numerimir  of  the  first  term  on  the  riglu-hand  skfe  of  Equation  (12)  is  die 
toj^gftadinal  foordl  central  mtHnent  of  the  concentration  distribution.  Mcnnents  of  cnder 
l^l^hbr  dam  fiwr  could  be  calculated,  but  were  not  needed  to  adequately  describe  the  plume 
b^vibr. 

Nuinriical  implementation  of  the  moments  calculation  was  necessarily  three 
^raenric^.  The  dcnnain  of  integration  was  discretized  into  volumetric  subdomains 
liadcSilBd  widi  the  indivkhial  sanqiling  poims.  A  primary  triangular  grid  was  constructed  in 
iSie  horizontal  {dane  using  the  MLS  as  ^id  points.  The  horizontal  region  of  each  subdomain 
by  constnicdng  a  secondary  mesh  with  vertices  at  the  centroids  of  triangles  and 
K  mid^ints  of  dieir  sides,  as  shown  on  Figure  17.  The  shaded  region  of  the  figure 
OMieiqponds  to  die  surface  area,  of  the  subdomain  associated  with  sampling  point  p,  and 
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F" 


is  e<pial  to  tbe  sum  of  one-third  of  the  areas  of  all  primaiy  grid  triangles  sharing  a  vertex  at  the 
point  p. 

The  vertical  dimension  of  the  subdomain,  Zp,  was  computed  as  the  sum  of  one-half  the 
distances  between  the  sampling  points  directly  above  and  below  p.  For  the  uppermost  and 
kiwmooat  sanqiling  points  on  each  MLS,  the  vertical  dimension  of  the  subdomain  was  set  equal 
to  cme-half  the  di^ance  between  p  and  the  next  vertically-adjacent  point,  plus  0. 19  meters  (i.e. , 
(Mie-half  the  nmmal  sample  point  spacing).  The  pore  fluid  within  each  subdomain  was  assunwid 
to  have  uniform  concentration  equal  to  the  measured  tracer  concentration  at  p.  The  mass  of 
tracer  associated  with  tbe  sample  point  p  is  defined  by: 

nip  =  nCpZpAp  (13) 

The  discrete  form  of  Equation  (5)  is, 
n 

Mijk=  EmpXp‘ypjZp^  (14) 

p-i 

where  n  is  die  total  number  of  sanqiliitg  points  and  Xp,  yp,  and  Zp  are  the  coordinates  of  point 
p.  The  accuracy  of  the  above  mass  integration  scheme  was  evaluated  through  mimepcal 
experimeids  in  which  calculations  were  compared  with  analytical  solution  results  for  a  test  plume 
c^.  Details  of  methods  and  results  of  this  evaluation  are  presented  in  Appendix  B. 

DATA  PREPARATION  FOR  MOMENTS  ANALYSIS 

In  the  few  instances  where  the  MLS  network  did  not  adequately  bound  the  plume 
laterally,  adjaceirt  imaginary  MLS  points  having  zero  concemration  were  added  at  6  meters  (the 
^kal  intrarow  MLS  placing)  outside  die  boundary  MLS  location  of  concern.  This  permitted 
limited  borizomai  extr^iolation  of  phime  measurements  outside  the  sampled  MLS  array. 
Votiad  extnqiolation  of  the  phime  was  also  performed  in  cases  where  the  uf^r  or  lower  point 
(m  an  MLS  stowed  above  background  (or  above  detection)  concentrations.  This  extrapolation 
was  based  on  imaginary  zero  concentration  points  added  at  0.38  meters  (one  typical  vertical 
sarr^  poirtt  placing)  above  os  below  the  non-zero  concentradon  boundary. 

Measui^  tritium  and  counts  were  corrected  for  background  acdvi^r.  Estimated  mean 
background  activities  of  2.0  and  2.8  pCi/mL  were  subtracted  from  the  tritium  and  '*C 
measurements,  reflectively.  The  tritium  data  were  corrected  for  radiological  decay,  but  the 
^*C  data  were  not  because  of  the  long  half  life  of  The  total  decay  during  the  440  day  field 
stiKly  was  about  6.7  percent  for  tritium  and  less  than  0.0002  percent  for  '*C. 

Two  corrections  were  applied  to  the  organics  and  data,  referred  to  as  Correction  I 
and  Cmnectimi  II.  Correction  I  was  ai^lied  to  data  sufiected  to  be  affected  by  cross¬ 
contamination  during  sampling.  Sample  data  exhibiting  relative  concentrations  of  '*C  or  an 
organic  compound  greater  than  the  relative  concentration  for  tritium  were  considered  suspect. 
In  general,  and  organic  compound  relative  concemrations  in  a  sample  should  be  less  than 
or  e^ial  to  the  cmre^nding  relative  concentration  for  tritium,  due  to  their  aisceptibility  to 
attenuating  processes  such  as  sorption,  volatilization,  and  biotransformation.  Suspect  ’^C  or 
mgamc  concentration  meamrements  were  adjusted  such  that  their  relative  concentrations  were 
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equal  to  that  of  the  corresponding  tritium  measurement  using  the  following  equation; 


C*i(org)  =  C„(org)Ci(’H)/C„(^H)  (15) 

In  this  equation,  C*i(org)  is  the  corrected  organic  or  concentration  for  sample  i,  Co(org)  is 
the  initial  injected  concentration  of  the  or  organic  compound,  CX^H)  is  the  tritium 
concentration  of  sample  i,  and  Co(^H)  is  the  initial  injected  tritium  concentration.  Evidence  of 
cross-contamination  was  limited  almost  entirely  to  data  from  the  first  snapshot  . 

Correction  n  was  aj^lied  to  samples  with  the  background  tritium  activity  level,  but 
having  detectable  amounts  of  one  or  more  of  the  organic  compounds.  For  these  samples,  the 
suspect  organic  solute  concentrations  were  set  to  zero.  This  procedure  eliminated  spurious 
analytical  results  in  the  low-concentration  range,  which  could  influence  the  higher  moments 
estimates. 

The  sensitivity  of  the  solute  mass  balance  estimates  to  these  corrections  is  illustrated  on 
Figure  18.  It  is  clear  that,  except  for  the  flrst  snapshot,  mass  estimates  for  the  four  organic 
solutes  are  insensitive  to  the  data  correction  method.  However,  this  is  not  the  case  for  ‘^C. 
After  132  days,  the  mass  estimates  of  were  highly  sensitive  to  application  of  these  data 
corrections,  due  to  the  presence  of  at  near-background  concentrations  in  most  samples.  On 
the  basis  of  these  results.  Correction  I  was  applied  to  the  organics  and  data  for  the  first 
snapshot.  Corrections  I  and  II  were  both  appli^  to  data  from  the  later  snapshots. 

SOLUTE  MASS  (ACTIVITY)  AT  SAMPLING  TIMES 

Relative  mass  (activity)  balances  for  each  solute  are  given  in  Tables  6-11,  along  with 
summaries  of  the  higher  moments  estimates.  Figure  19  presents  the  relative  mass  balance  trends 
during  the  experiment.  No  relative  activity  or  higher  moment  estimates  are  given  for  tritium 
or  '^C  for  the  snapshot  at  440  days.  As  noted  earlier,  sampling  at  440  days  was  limited  to  the 
near-field  region,  and  thus  did  not  completely  encompass  the  radiological  tracer  plumes  and  may 
have  only  partially  bounded  the  o-DCB  plume. 

The  activity  balance  for  tritium  shows  an  overestimate  at  27  days,  followed  by 
approximately  100  percent  recovery  for  the  two  succeeding  snapshots,  and  77  percent  recovery 
for  the  fmal  smq>shot  at  328  days.  The  initial  overestimation  of  tritium  activity  is  attributed  to 
the  same  sampling  artifacts  that  produced  the  excess  bromide  recovery  during  the  early  stages 
of  MADE-1 .  That  is,  preferential  sampling  by  the  MLS  from  relatively  permeable  aquifer  zones 
in  which  solute  concentrations  are  locally  high  may  lead  to  mass  overestimation  because  the 
mass  integration  scheme  assumes  local  uniformity  of  tracer  concentration  (Boggs  and  Adams, 
1992).  In  addition,  suspected  vertical  hydraulic  interconnection  among  sampling  points  associated 
with  certain  MLS  located  close  to  the  injection  site  may  result  in  over-pr^iction  of  tracer  mass 
(Boggs  et  al.,  1988).  The  23  percent  tritium  deficit  at  328  days  was  due,  at  least  in  part,  to 
apparent  frontal  truncation  of  plume  which  had  then  passed  the  outermost  boundary  of  the 
well  fleld. 

Aside  from  the  frontal  truncation,  lateral  bounding  of  the  remainder  of  the  tritium  plume 
appears  to  have  been  reasonably  complete.  As  evidence,  the  addition  of  the  imaginary  zero- 
concentration  MLS  adjacent  to  boundary  MLS  exhibiting  above-background  tritium  levels 
produced  an  increase  in  relative  activity  of  2  percent  or  less  for  each  snapshot.  Some  degree 
of  vertical  truncation  of  the  plume  was  indicated  by  the  occurrence  of  above-background  tritium 
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Legend 

— ^  uncorrected 
correction  I 
»  correction  Q 


Figure  18.  Sensitivity  of  Organic  Solute  Mass  Estimates  to  Data  Preparation  Method. 
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TABLE  6.  TRITIUM  PLUME  CHARACTERISTICS 


Tiin^  (<lays) 


27 

132 

224 

328 

440t 

M/M„ 

1.52 

1.05 

0.98 

0.77 

— 

Xc  (m) 

0.0 

-0.9 

0.2 

2.1 

— 

yc  (m) 

3.9 

8.1 

46.5 

76.8 

— 

zc  (m) 

58.22 

58.68 

58.80 

58.48 

— 

Hor.  Displ.  (m) 

3.9 

8.2 

46.5 

76.8 

— 

Ver.  Displ.  (m) 

0.42 

0.88 

1.00 

0.68 

— 

On  (m^) 

10.3 

94.4 

4,380 

6,560 

— 

o-n  (m^) 

8.6 

7.9 

52.5 

107. 

— 

Oi^  (m^) 

2.0 

1.3 

2.5 

2.9 

— 

gi 

0.72 

1.5 

0.74 

0.32 

— 

g2 

1.5 

6.3 

0.34 

-0.58 

— 

©h  (deg.  CCW  from  y-axis) 

11.1 

0.2 

2.3 

3.3 

— 

(deg.  CCW  from  y-axis) 

-1.1 

3.1 

1.7 

2.0 

_ 

t  moments  not  estimated  for  spatially  limited  sampling  at  440  days 
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TABLE  7.  BENZENE  PLUME  CHARACTERISTICS 

Tiny  fflayf?) 


27 

132 

224 

328 

440 

M/Mo 

0.92 

0.43 

0.23 

0.07 

0.06 

Xc  (m) 

-0.2 

-1.1 

-0.8 

-1.0 

-1.1 

yc  (m) 

3.8 

6.3 

12.4 

7.7 

7.9 

Zc  (m) 

58.13 

58.68 

59.24 

58.90 

58.67 

Hor.  Displ.  (m) 

3.8 

6.4 

12.4 

7.7 

7.9 

Ver.  Displ.  (m) 

0.33 

0.88 

1.44 

1.10 

0.87 

(m^) 

9.2 

38.4 

826 

24.7 

21.4 

^22*  (m^) 

7.9 

6.9 

6.1 

8.4 

10.7 

<^33^  (m^) 

1.9 

1.1 

1.2 

1.2 

1.3 

gi 

0.67 

2.6 

2.2 

0.61 

0.08 

gz 

1.4 

25.8 

8.9 

1.6 

-0.13 

l^h  (<leg.  CCW  from  y-axis) 

39.0 

0.4 

3.6 

1.6 

14.1 

6y  (d^.  CCW  from  y-axis) 

3.8 

6.7 

-3.7 

1.7 

-4.4 
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TABLES.  NAPHTHALENE  PLUME  CHARACTERISTICS 


liffw  (tfays) 


27 

132 

224 

328 

440 

M/M, 

0.58 

0.45 

0.25 

0.08 

0.06 

Xc  (m) 

-0.3 

-1.2 

-1.4 

-1.0 

-1.8 

yc  (m) 

3.5 

5.8 

6.6 

7.2 

7.3 

Zc  (m) 

58.05 

58.54 

59.03 

58.71 

58.62 

Hor.  Displ.  (m) 

3.5 

5.9 

6.8 

7.2 

7.6 

Ver.  Displ.  (m) 

0.25 

0.74 

1.23 

0.91 

0.82 

<^11^  (m^) 

8.5 

19.2 

14.7 

16.2 

12.8 

022^  (m^) 

7.3 

6.7 

6.1 

4.8 

4.7 

(m^) 

1.6 

1.4 

1.3 

1.0 

1.3 

gi 

0.77 

0.52 

0.44 

0.17 

0.15 

g2 

1.7 

0.24 

0.29 

-0.43 

-0.41 

(deg.  CCW  from  y-axis) 

84.5 

-9.8 

-5.5 

0.5 

5.5 

$y  (deg.  CCW  from  y-axis) 

-0.8 

-2.0 

-1.4 

1.1 

-16.5 
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TABLE  9.  p>XYLENE  PLUME  CHARACTERISTICS 


Tfang  fdayrt 


27 

132 

224 

328 

440 

M/Mo 

0.77 

0.58 

0.23 

0.03 

0.01 

Xc  (m) 

-0.2 

-1.0 

-1.2 

-1.1 

-1.3 

yc(m) 

3.6 

6.1 

6.4 

10.5 

6.4 

Zc(m) 

58.04 

58.60 

58.93 

58.66 

57.93 

Hor.  DU^l.  (m) 

3.7 

6.1 

6.5 

10.5 

6.6 

Ver.  Di^l.  (m) 

0.24 

0.80 

1.13 

0.86 

0.13 

ou  (m^) 

9.1 

22.8 

16.3 

213. 

11.8 

<^22^  (m^) 

7.2 

5.8 

6.4 

4.6 

4.0 

ajj^  (m^) 

1.9 

1.1 

1.3 

1.0 

1.3 

gi 

0.82 

0.79 

0.81 

1.7 

-0.12 

g2 

1.9 

3.1 

3.3 

5.2 

-0.72 

(d^.  CCW  ftcm  y'ttis) 

13.2 

-11.3 

5.0 

3.4 

28.0 

dv  (deg.  CCW  from  y-axis) 

5.3 

-1.8 

6.6 

1.9 

2.2 
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TABLE  10.  o-DCB  PLUME  CHARACTERISTICS 


Time  (davs) 


27 

132 

224 

328 

440 

M/Mo 

0.81 

0.80 

0.60 

0.21 

0.13 

Xc  (m) 

0.0 

-1.0 

0.0 

-0.2 

-1.5 

yc(m) 

3.9 

7.4 

34.7 

22.1 

8.1 

Zc  (m) 

58.09 

58.57 

58.93 

58.90 

58.68 

Hor.  Displ.  (m) 

3.9 

7.5 

34.7 

22.1 

8.2 

Ver.  Displ.  (m) 

0.29 

0.77 

1.13 

1.10 

0.88 

ou  (m^) 

11.3 

64.4 

3540 

1180 

21.9 

On  (m^) 

7.9 

7.8 

25.2 

22.4 

7.1 

On  (m^) 

2.1 

1.4 

2.0 

1.6 

1.8 

8i 

0.82 

1.7 

1.0 

1.6 

0.12 

82 

1.6 

8.3 

1.2 

6.5 

-0.22 

Ofc  (deg,  CCW  from  y-axis) 

4.9 

0.5 

2.5 

5.5 

7.5 

(?y  (deg.  CCW  from  y-axis) 

-1.6 

5.5 

2.4 

0.7 

-5.4 
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tASLBll.  CARBON-14  PLUME  CHAIUCTERISTICS 

Thwfe  rdavw^ 


27 

132 

224 

328 

440t 

M/H> 

0.72 

0.65 

0.49 

0.40 

Xc  (m) 

-0.2 

-1.0 

-0.4 

1.5 

“ 

yc(*n) 

4.1 

6.7 

30. 

68.5 

— 

Zc  (n») 

58.22 

58.60 

58.79 

58.47 

“ 

Hor.  Displ.  (m> 

4.1 

6.8 

30.0 

68.6 

— 

Ver.  D^l.  (m> 

0.42 

0.80 

0.99 

0.67 

“ 

<t\\  (n^) 

11.4 

39.8 

2930 

6490 

— 

(m^ 

8.6 

7.5 

34.7 

97.6 

ff33*(in*) 

1.9 

1.4 

2.3 

2.7 

“ 

gi 

0.71 

1.4 

1.2 

0.42 

— 

1.3 

6.2 

2.3 

-0.46 

— 

^  (ddg.  CCW  fiom  y-axis) 

9.1 

1.3 

1.8 

2.9 

~ 

(dq{.  CCW  from  y-axis) 

-0.8 

6.3 

0.8 

2.1 

_ 

t  moments  not  estimated  for  ^tially  limited  sampling  at  440  days 
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Figure  19.  Benzene,  Naphthalene,  p-Xylene,  and  o-DCB  Relative  Mass  Balances. 
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BMMMWMeKi  at  die  uppermosc  and  lowermost  san^ling  poims  of  certain  MLS.  The  precise 
inagailMte  of  (he  bmos  (activity)  recovery  error  involved  in  die  vertical  tnincatkm  is  unknown. 
Hoaeevor,  aendlivity  raises  in  whidi  5ctitious  zero^oncentntfion  points  were  added  at  a 
diiMKeof  1.14  aaeim  (three  sample  point  spacings)  instead  of  the  normal  0.38  meters,  indicate 
amis  recovery  increases  nmgmg  ftom  7  pacen.  of  injected  mass  for  snqidiot  1  to  2  percent  for 
anapahot  4.  These  residts  mate  it  unlikely  that  die  magnitude  of  the  errm’  associated  widi  die 
mna  recovery  eathnales  is  greater  than  tewt  10  percent  of  the  injected  mass.  This  conclusion 
tt  consistent  wift  ttmt  thawn  by  Adams  and  Gelhar  (1992)  regarding  vertical  truncation  errors 
associated  with  tte  MADE-1  Inomide  experimmit. 

The  otgaiuc  sohnes  exhibit  decreasing  mass  recoveries  throughout  die  field  experiment. 
Dataaeiie  shows  an  abrupt  decline  in  mass  recovery  between  27  and  132  days,  followed  by  a 
deoeaidng  rate  of  mass  decline  during  the  renuiinider  of  the  study.  Excqit  for  a  more  gtadiuil 
mhiil  BMM  decrease,  the  nsass  balance  trends  for  nqihdudene  and  p-xylene  ate  similar  to  that 
of  benzene.  Hie  relative  mass  of  o-DCTB  is  essentially  stable  duriiig  die  first  132  days  after 
whidi  an  mcreasing  rate  of  nms  loss  is  observed.  The  appareitt  order  of  organics  removal  was, 
ftom  gremest  to  least,  benzene  fdlowed  by  p-xylene,  naidtthaloK,  and  o-DCB.  After  440  days, 
die  rdadve  mass  recoveries  for  benzene,  naididialene,  p-xylene,  and  o-DCB  were  6,  6,  1,  and 
13  percent,  respectively.  The  fmal  mam  estimate  for  o-DCB  probably  represeitts  a  lowcx'  limit, 
smce  some  &egttt  of  dimmtkMi  of  die  o-DCS  plume  was  evident  in  the  sample  measurements. 

The  dispathy  between  the  mass  balan^  for  p-xylene  and  is  noteworthy.  The  relative 
iMweit  Of  p-xylene  and  were  tpuie  similar  at  27  days  (i.e.,  77  and  72  percent,  respectively). 
The  lidtial  simiarii^  is  expected  since  was  chmically  associated  with  p-xylene  at  die 
hqginiihig  xif  the  expetriment.  However,  the  rviladve  mass  curves  for  the  two  compounds 
sttbaeguent^  ^voge  mdiciuing  transfmniuxni  of  die  p-xylene.  After  328  days,  the  relative 
mass  recovered  ftR*  is  40  percent,  vdiile  diat  of  p-xylene  is  only  3  petcem. 

CENTER  CF  MASS  AND  VELOCITY 

Fipiie  20  shows  die  raemi  horizmnal  trajectories  of  the  solute  plumes  as  indicated  by  die 
norttml  projections  cf  die  first  moments  of  die  solute  plumes  on  a  horizontal  plane.  Duriiig  the 
expohnent,  die  tritium  phune  followed  an  approxiimuely  linear  trajectory,  which  was  nearly 
puiHel  to  positive  y-axis  of  die  site  coordinate  system.  This  overall  direction  of  movement 
b  cottsistmit  with  the  horizcmtal  hydraulic  gradient  shown  earlier  in  Figure  3.  exhibited  a 
najectory  similar  to  that  of  tritium,  but  with  a  slighdy  lower  overall  displacement.  The 
hcMtkcnital  trajectories  of  the  other  solutes  were  similar  to  that  of  tridum  durmg  the  first  132 
days  of  die  study,  but  deviated  markedly  diereaft^.  All  of  the  organic  ctunpounds  had  smaller 
mmn^liS|daoemeilts  dum  did  tritium,  and  an  anwrent  "retreat"  of  all  organic  plumes  except 
B  indtoded  during  (he  latter  stages  of  the  experiment.  This  retreat  is  a 
tnaaifbStation  of  mgiinic  solute  degradation. 

Figure  21  shows  the  mean  horizontal  displacement  of  the  solute  plumes  during  the  study 
period.  The  mean  velocity  indicated  for  tritium  durii^  the  first  132  days  is  23  m/yr,  while  the 
mom  velocity  during  the  ranainder  of  the  study  is  128  m/yr.  The  increase  in  the  mean  velocity 
is  conskitent  widi  the  increasing  mean  hydraulic  ctmductivity  along  the  phime  travel  path,  and 
reflects  the  increasing  fraction  of  tritium  entering  the  relatively  high-velocity  gnxmdwater 
present  in  the  far-fiekl.  Naphthalene  displacement  increased  monotonically  over  time,  but  at  an 
iqiparent  rate  of  oidy  6.8  ni/yr.  The  mean  horizontal  displacements  of  the  benzene,  p-xylene. 
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Figure  20.  Solute  Plume  Trajectories  Shown  as  Normal  Projections  of  Center  of  Mass  Locations 


Days 


F^vie  21.  Tntqmral  Variation  in  (a)  Relative  Mass  Balance,  (b)  Mean  Horizontal 
Di^lacement  and  (c)  Mean  Vertical  Displacement,  for  Each  Solute. 
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and  o-DCB  plumes  decreased  during  the  latter  stages  of  the  study.  The  displacement  of 
was  slightly  less  than  that  of  tritium,  but  much  greater  than  that  of  p-xylene. 

To  better  quantify  the  average  horizontal  groundwater  velocity  of  the  near-field, 
where  the  majority  of  tritium  tracer  resided  during  the  experiment,  a  separate  set  of  Hrst 
moments  was  estimated  using  only  the  tritium  measurements  for  MLS  located  within 
ai^roximately  25  meters  of  the  injection  point.  These  results  are  presented  on  Figure  22. 

The  mean  velocity  indicated  from  the  tritium  displacement  between  27  days  and  440  days  is 
approximately  5  m/yr  (0.014  m/d).  The  vertical  positions  of  the  tracer  plume  centroids 
during  the  experiment  are  shown  on  Figure  21.  All  plumes  exhibited  upward  movement 
during  the  first  224  days,  followed  by  downward  or  near-stable  vertical  displacement.  The 
initial  upward  displacement  reflects  the  upward  vertical  hydraulic  gradients  present  in  the 
near-field.  The  relatively  low  centroid  elevations  indicated  for  tritium  and  after  224  days 
are  attributed  to  the  presence  of  a  substantial  fraction  of  these  solutes  in  the  far-field,  where 
hydraulic  gradients  were  downward  to  near  neutral.  No  evidence  of  vertical  movement 
resulting  ^m  densiQr  differences  between  the  injection  solution  and  ambient  groundwater  is 
apparent  in  these  results.  A  density-induced  hydraulic  gradient  of  O.OOOIS  cm/s  in  the 
downward  direction  is  estimated  for  the  tracer  injection  solution.  The  density  gradient  is 
small,  conq>ared  to  the  mean  upward  vertical  gradient  of  0.04  cm/s  measured  near  the 
injection  site. 

SPATIAL  VARIANCE,  SKEWNESS,  AND  KURTOSIS 

Figure  23  shows  the  principal  components  of  the  spatial  variance  of  the  tritium  plume 
as  a  function  of  mean  horizontal  displacement.  The  longitudinal  variance  trend  is 
approximately  linear  after  the  first  10  meters  of  displacement.  The  slight  reduction  in  the 
rate  of  increase  indicated  between  displacements  of  46.S  and  76.8  meters  is  attributed  to 
frontal  plume  truncation.  The  horizontal  and  vertical  transverse  variances  both  exhibit  slight 
decreases  in  magnitude  between  3.9  and  8.2  meters,  followed  by  gradual  nonlinear  increases 
out  to  a  di^lacement  of  76.8  meters. 

Temporal  trends  in  the  longitudinal  variance,  skewness,  and  kurtosis  for  all  solutes 
are  compared  on  Figure  24.  The  higher  moments  are  useful  to  the  extent  that  they  further 
emphasize  some  of  the  similarities  and  differences  in  the  ensemble  characteristics  of  the 
solute  plumes.  Similar  trends  for  the  higher  moments  are  observed  for  and  tritium.  Both 
exhibit  nonlinear  increases  in  longitudinal  variance,  and  both  show  increasing  skewness  and 
kurtosis  up  to  132  days,  followed  by  a  gradual  decrease  to  near  zero  at  440  days.  The  four 
organic  solutes  exhibit  widely  vaiying  higher  moment  trends  apparently  due  to  differential 
transformation  and  sorption  effects.  The  higher  moments  associated  with  the  o-DCB  plume 
are  similar  to  those  for  tritium  up  to  224  days,  but  they  then  deviate  subsequently  from 
tritium  due  to  the  diminishing  longitudinal  extent  of  the  o-DCB  plume.  The  higher  moments 
for  benzene,  naphthalene,  and  p-xylene  are  similar  to  those  of  tritium  only  for  the  snapshot  at 
27  days,  after  which  they  are  quite  dissimilar  from  tritium  and  from  each  other. 
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Days 


Figure  22.  Temporal  Variation  of  Mean  Horizontal  Displacement  of  the  Tritium  Plume 
Based  on  Spatial  Moments  Analysis  of  Near-Field  Data. 
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Figure  24.  Temporal  Trends  in  (a)  Longitudinal  Spatial  Variance 
(b)  Skewness  and  (c)  Kurtosis,  for  Each  Solute. 


SECTION  Vi 

HYDROCARBON  DEGRADATION  RATES 


Two  separate  sets  of  laboratory  studies  were  performed  prior  to  the  field  experiment 
in  order  to  predict  the  degree  of  degradation  of  the  organic  compounds  that  might  occur  in 
the  field.  It  was  understood  from  the  inception  of  these  studies  that  laboratory  measured 
degradation  rates  would  probably  be  different  than  field  degradation  rates,  and  that  laboratory 
results  could  identify,  but  not  quantify  the  role  of  degradation  in  the  field  experiment. 
Procedures  and  results  of  the  laboratory  studies  are  given  below. 

MIXED  BATCH  REACTOR  STUDIES  OF  DEGRADATION 

Mixed  batch  reactor  transformation  studies  were  conducted  with  Columbus  aquifer 
material  to  (1)  determine  the  potential  for  degradation  of  the  three  aromatic  hydrocarbon 
compounds  selected  for  the  field  experiment  (i.e.,  benzene,  naphthalene,  and  p-xylene),  and 
(2)  to  qualitatively  compare  the  degradabiliQr  of  the  more  biologically  resistant  o-DCB  with 
that  of  the  aromatic  hydrocarbons.  These  studies  were  qualitative  in  nature,  and  were  not 
intended  to  provide  quantitative  biotransformation  data  for  comparison  with  the  field 
experiment. 

Approximately  20,000  kg  of  aquifer  material  was  excavated  from  the  saturated  zone 
of  the  Columbus  aquifer  to  facilitate  these,  and  other  laboratory  tests.  The  excavation 
location  was  just  outside  of  the  well  field  on  the  southwest  comer  of  the  test  site.  A 
subsample  of  this  material  was  wet  sieved  on-site  with  groundwater  using  a  2  mm  sieve  and 
transferred  to  a  cooler.  To  minimize  biological  contamination,  sieves  and  coolers  were 
cleaned  with  ethanol  and  air-dried  prior  to  use.  The  samples  were  then  saturated  with 
groundwater,  packed  in  ice,  and  transported  to  the  AL/EQC  Laboratories  at  Tyndall  AFB 
where  they  were  stored  at  S^C.  Approximately  40  liters  of  groundwater  was  also  collected 
for  use  in  these  studies. 

Saturated  aqueous  solutions  of  benzene,  p-xylene,  and  naphthalene  in  Milli-Q  water 
were  prepared  in  sterilized  saturator  flasks  as  described  by  Burris  and  MacIntyre  (1985). 
Aliquots  of  these  saturated  solutions  were  diluted  with  site  groundwater  to  produce 
hydrocarbon  concentrations  of  approximately  10  mg/L.  These  solutions  were  added  to 
sterilized,  10  mL  vials  containing  approximately  10  grams  of  aquifer  material.  Active  vials 
were  then  sealed  with  Teflon*-backed  septums  and  aluminum  crimp  seals.  Poisoned  vials 
were  spiked  with  sodium  azide  to  bring  the  azide  concentration  to  0.02  percent  to  inhibit 
degradation,  and  they  were  also  sealed.  Finally,  control  samples  containing  azide  and 
hydrocarbons,  but  no  aquifer  material,  were  placed  in  vials  and  sealed.  Each  of  the  three 
different  treatments  vvas  repeated  in  triplicate.  The  sealed  bottles  were  placed  on  a  rotator 
spinning  at  30  revolutions  per  minute.  The  rotator  was  located  inside  an  incubator  held  at 
25 “C.  Triplicate  samples  of  each  treatment  were  removed  from  the  rotator  at  weekly 
intervals.  Sample  bottles  were  centrifuged  at  2,000  rpm  and  3.0  mL  of  aqueous  hydrocarbon 
solution  was  removed  and  added  to  a  4.0  mL  vial.  One-half  mL  of  pentane,  containing 
toluene  as  an  internal  standard,  was  then  added  to  the  vial  and  the  hydrocarbons  were 
extracted  from  the  aqueous  phase. 
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S^^nration  and  analysis  of  the  extracted  hydrocarbons  was  carried  out  on  a  Hewlett- 
Packard  Model  85  gas  chromatograph  (GC)  equipped  with  a  30-meter  c^illaiy  colunm 
coated  with  methyl  sUicone.  The  C^C  was  temperature  programmed  for  an  initial  2-minute 
hold  at  40*’C,  followed  by  a  steady  increase  to  ISO'^C  at  lS°C/min,  which  was  held  for  6 
minites.  Flame  ionization  detection  was  used  to  quantify  the  separated  hydrocarbons. 
Relative  refuse  factors  of  the  hydrocarbons  were  determined  so  that  final  concemrations 
could  be  calculated  from  the  peak  areas  of  the  separate  samples  and  the  toluene  internal 
standard. 

Results  from  this  set  of  experiments  (designated  as  Study  1)  showed  that  the 
hydrocarbons  disappeared  in  the  active  vials  (Figure  25),  but  remained  relatively  constant  in 
fbc  vials  containiqg  sodium  azide  inhibitor.  Control  solutions  showed  a  slow  loss  over  the 
eiqierimeittal  period.  These  losses  may  be  attributed  to  sorption  and/or  vaporization,  but  are 
small,  conqiated  to  the  rapid  loss  of  hydrocarbons  in  the  active  vials. 

Due  to  the  complete  disappearance  of  hydrocarbons  in  Study  1,  a  second  set  of 
experiments  was  conducted.  Study  2  was  a  complete  replication  of  Study  1,  except  that 
o-DCB  was  added  to  the  suite  of  compounds.  This  chemical  is  a  common  industrial  solvent 
ftHind  at  numerous  groundwater  contamination  sites,  and  was  purposely  selected  for  its 
greater  resistance  to  ttegiadation. 

Results  from  Study  2  (Figure  26)  are  similar  to  those  of  Study  1  except  that  the 
o-DCB  showed  a  greater  resistance  to  degrsulation  than  the  other  test  compounds,  as  was 
anticq)ated.  There  were  measurable  quantities  of  o-DCB  left  in  the  active  vials  after  15 
days.  Based  on  the  results  from  these  preliminary  experiments,  it  was  decided  to  include 
o-DCB  in  the  suite  of  conqwunds  to  be  used  in  the  field  test,  to  provide  an  organic 
conqxxind  with  a  potentially  lower  degradation  rate. 

The  results  of  these  experintents  were  useful  in  designing  the  field  experiment,  and 
indicated  that  degradation  of  the  selected  aromatic  hydrocarbons  in  the  aquifer  might  be  quite 
nqiid.  The  investigators  proceeded  with  the  fkld  experiment  under  the  assumption  that 
degradation  of  these  sohnes  would  occur  in  the  aquifer,  but  at  a  considerably  slower  rate 
than  in  the  laboratory.  Fortunately,  this  assumption  proved  to  be  a  correct  one. 

DEGRADATION  OF  ORGANICS  FROM  MOMENTS  ANALYSIS 

As  described  in  Section  4,  dissolved  oxygen  samples  were  taken  8  days  before  and 
48,  111,  161,  264  and  330  days  after  the  test  injection.  These  samples  were  taken  at  three 
discrete  levels  from  a  set  of  16  multilevel  sampler  wells  located  along  the  plume  axis, 
including  wells  at  2.3  meters  upgradient  and  3.6  meters  down  gradient  of  the  injection  well?. 
The  results  show  that  the  organic  solutes  degraded  under  aerobic  conditions  in  the  aquifer. 
This  outcome  was  anticipated  because  the  initial  organic  solute  concentrations  in  the  aquifer 
after  mixing  were  estimated  to  be  too  low  to  signiHcantly  deplete  the  existing  level  of 
dissolved  oxygen  near  the  injection  wells.  These  DO  concentrations  were  high  enough  so 
that  degradation  rates  of  the  organic  compounds  were  probably  not  oxygen-limited.  In 
support  of  this  assumption,  Hlinel  (1988)  notes  that  the  dissolved  oxygen  content  in  a 
ncmnally  (grating  activated  sludge  plant  is  1-2  mg/L,  which  is  adequate  to  support  aerobic 
biodegradation  processes. 

Figure  2  shows  the  distribution  over  a  vertical  section  oriented  along  the  test 
phime  axis,  and  indicates  large-scale  heterogeneity  and  structures  that  Rehfeldt  et  al.  (1992) 
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Figure  25.  Batch  Reactor  Results  from  Study  1. 
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Fi^te  26.  Batch  Reactor  Results  from  Study  2. 
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refer  to  as  channels.  Close  to  the  injection  wells  is  relatively  low,  with  values  near  10'^ 
cm/s  in  a  region  extending  from  the  phreatic  surface  to  the  lower  confining  layer,  and  to 
ai^roximately  40  meters  down-gradient  from  the  injection  location.  From  40  meters  out  to 
200  meters,  K,.  in  the  upper  3  meters  of  the  aquifer  is  about  10'‘  cm/s,  while  in  the  deeper 
portion  of  the  aquifer  remains  low.  Thus,  a  near  surface  channel  crosses  the  intermediate 
and  far  fields. 

Most  of  the  injected  solute  mass  remained  in  the  low  region  (near  the  injection 
wells)  throughout  the  course  of  the  experiment.  This  situation  allowed  degradation  rate 
constants  to  be  determined  separately  from  data  for  the  entire  plume,  and  from  data  for  only 
the  near  field.  Solute  reactions  included  tritiated  water  (Tr)  decay,  with  a  of  1.548  x  10^ 
days  ’,  and  degradation  of  organic  compounds  with  rate  constants  k^,  k^,  k„  and  k^  for 
benzene,  p-xylene,  naphthalene  and  o-dichlorobenzene,  respectively.  Degradation  rate 
constants  based  on  data  from  the  whole  field  were  calculate  for  each  snapshot  using  the 
zero"*  moments  of  the  plumes  of  tritiated  water  and  of  the  organic  solutes.  For  each  organic 
solute,  initial  mass  and  mass  remaining  for  each  snapshot  sampling  time  given  in  Tables  7-10 
were  used  in  a  nonlinear  least-squares  regression  to  fit  a  first-order  exponential  function,  thus 
assuming  a  flrst-order  degradation  process.  Degradation  rate  constants  from  these 
regressions  are  given  in  Table  12.  Maximum  degradation  rates  for  all  solutes  occurred 
between  224  and  328  days  after  injection,  as  can  be  seen  from  Figure  27,  and  these  rates  are 
included  as  k„„,  in  Table  12. 


TABLE  12.  DEGRADATION  RATE  CONSTANTS  FOR  ORGANIC 
COMPOUNDS  USING  WHOLE-FIELD  DATA 


k  (days  ’) 

lima,  (days  ’) 

benzene 

0.0070 

0.011 

p-xylene 

0.0107 

0.020 

naphthalene 

0.0064 

0.011 

o-dichlorobenzene 

0.0046 

0.010 

The  observed  mass  loss  is  presumed  to  be  due  to  microbial  catabolism 
(biodegradation),  but  is  referred  to  simply  as  degradation  since  biodegradation  is  indicated, 
but  not  proven  by  this  experiment.  Proof  of  biodegradation  in  the  Columbus  aquifer  would 
require  field  application  of  molecular  biology  techniques  similar  to  those  described  by 
Fleming  et  al.,  (1993). 

Spatial  distributions  of  solutes  are  apparently  related  to  the  distribution,  as 
indicated  by  Figure  28(a),  which  shows  the  distribution  of  relative  tritium  concentrations  at 
224  days  after  injection.  High  concentrations  remained  near  tlie  injection  wells,  with 
lower  concentrations  in  the  upper  portion  of  the  aquifer  from  about  30-250  meters 
downgradient.  There  was  little  in  the  lower  portion  of  the  aquifer  between  30  and  175 
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DAYS  AFTER  INJECTION 


Figure  27.  Plot  of  Organic  Compound  Masses  Versus  Time, 
Based  on  Data  from  the  Whole  Field. 
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Figure  28.  Distribution  of  (a)  Relative  Tritium  Concentration  and  (b)  Relative  p-Xylene 

Concentration  Over  a  Vertical  Section  Along  the  Plume  Centerline  at  224  Days. 
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meien  downpadient.  This  distribution  suggests  that  most  of  the  tritium  was  transited 
from  die  low  zone  near  the  injection  wells  to  the  far  field  through  an  upper  chimnel  whose 
size,  locatkm,  and  govern  the  solute  loss  rate  from  the  near  field  region. 

The  organic  sohnes  were  slightly  retained,  and  thus  more  confined  to  die  near  field 
din  was  die  umetained  tritium.  Accordingly,  most  degradatimi  of  the  organic  solutes 
occurred  in  the  near  field  zone.  This  zone  is  qifHoximately  delineated  by  the  p-xylene 
distribidion  at  224  days  after  injection,  shown  in  Figure  28(b).  Most  of  the  p-xylene 
remained  within  20  meters  of  die  iiyection  wells,  and  the  othor  organic  solutes  were  similarly 
distributed.  They  were  primarily  confined  to  the  near  field  zone  because  they  degraded 
raindly  relative  m  their  transport  from  the  near  field.  As  a  result,  separate  degradation  rates 
could  be  cateulated  using  only  the  near  field  data,  treating  the  near  field  region  zone  as 
mixed  reactor  zone  with  a  slow  leak  to  the  far  field.  A  schematic  diagram  of  the  transport 
and  degradation  processes  occurring  in  the  near  field  region  is  shown  in  Figure  29.  Separate 
d^radation  rams  calculated  from  the  reactor  model  were  useftil  because  the  solute 
ccmcentratioiis  were  much  higher  in  the  near  field  and  less  subject  to  analytical  uncertainty. 

Sohtte  mass  transit  from  the  near  field  to  the  far  field  is  assumed  to  be  a  first-order 
process  with  the  same  constant,  k|,  for  all  solutes.  Degradation  and  loss  to  the  far  field 
occurred  simultaneously  in  the  near  field.  The  rate  etpiations  are; 


-d[Tr]/dt  =  kd[Tr]  +  k,ITr] 

(16) 

-dlBj/dt  =  k,tB)  -»■  k,{Bl 

(17) 

-d[Xl/dt  «  k,IX]  +  k,IX) 

(18) 

-dlNJ/dt  =  k.IN]  k,fNl 

(19) 

-d(q/dt  =  k,[C]  +  k,[CJ 

(20) 

The  initial  ccmditions  for  these  equations  are;  =  0;  [Tr]o  =  0.539  Ci;  [B]o  =  660  g; 
Pqo=  402g;  [N]o  =  70g;  (Qo  =  318  g. 

Integration  of  rate  Equations  (16)  to  (20)  with  these  conditions  gives  solute  masses  in 
the  near  field  as  a  function  of  time.  The  transport  rate  constant,  k,  is  calculated  from  the 
integrated  form  of  Equation  (16),  using  the  known  k^  of  tritium.  Degradation  constants  for 
die  organic  solutes  are  then  (tetermined  by  substitution  of  k,,  and  solving  the  integrated  forms 
of  Equations  (17)  to  (20)  for  k^,  k,,  k,,  and  kc,  respectively. 

Sorption  of  organic  solutes  was  not  included  in  rate  calculations  using  the  near  field 
^ta.  Aromatic  solutes  were  weakly  sorbed  on  Columbus  aquifer  material  (Lion  et  ai,  1990; 
MacIntyre  et  oL,  1991).  Using  subsamples  of  a  20,000  kg  composite  sample  of  Columbus 
aquifer  material,  batch  sorption  coefficients  measured  for  naphthene,  o-DCB,  p-xylene  and 
benzene  were  0.(^5,  0.065,  0.048  and  0.059  L/Kg,  respectively.  Retardation  effects  on 
orgaiBC  sohite  distributiom  in  the  near  field  were  not  detected  due  to  the  predominant 
infioence  of  degradation.  Sorption-desorption  processes  at  aquifer  material  surfaces  are 
aiMimiMi  to  be  fast  relative  to  the  chemical  degradation  processes  considered  here,  as 
indicaied  by  Spiro  (1989),  and  to  be  rate-limited  by  physical  transport  of  organic  solute 
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Figure  29.  Schematic  Diagram  of  Transport  and  Degradation  Processes  in  the  Near-Field  Region. 


nM^ecales  to  and  fttMn  the  solid-solution  interface.  As  was  noted  by  Robinson  et  aL  (1990) 
tor  tt^iene  d^radation  on  surface  soils,  the  aromatic  compounds  may  be  microbially 
degtaded  from  solution,  sorbent  surfaces,  or  both.  This  question  is  not  resolved,  but 
rotiNkm  processes  probably  did  not  limit  microbial  degradation  of  organic  solutes  due  to  the 
small  80ft)tion  coefficients  measured  in  this  experiment. 

of  organic  solutes  by  transport  (hnn  the  near  to  the  far  field,  was  assumed  to  be 
proportimial  to  the  loss  of  tritium.  The  near  field  is  defined  as  the  portion  of  die  plume 
volume  within  10  meters  dowqgradient  from  the  injection  wells.  The  amount  of  tritium  in 
die  near  5eld  at  a  given  time,  was  calculated  by  the  ^tial  integration  method  d^ribed 
above.  The  level  of  ^HjO,  corrected  for  decay,  decreased  exponentially  with  time,  as  shown 
in  Figure  30.  This  plot  indicates  that  solute  leakage  from  the  reactor  was  apparenUy  first 
order,  widi  leakage  rate  constant,  k,  =  S.4S  x  10^  days  '.  A  mass  balance  for  tritium  over 
the  entire  plume  showed  that  the  injected  tritium  mass  was  accounted  for  at  the  times  given 
in  Figures  27  and  31.  This  result  implies  that  the  tritium  mass  loss  in  the  near  field  was 
indeed  due  to  leakage. 

Degradation  rates  of  organic  solutes  in  the  near  field  are  shown  in  Figure  31 .  The 
curves  shown  in  Figures  27  and  31  suggest  that  degradation  in  the  Columbus  aquifer  may  be 
first  order,  with  an  initial  lag  period.  The  departure  of  the  curves  in  Figures  27  and  31  from 
lineaHty  is  attributed  to  microbial  adaimition,  cell  growth  and  substrate  limitation.  Spain 
(1990)  has  described  lag  period  dependence  on  microbial  growth  and  ad^itation  m  die 
substrate.  The  observed  lag  period  is  followed  by  degradation  at  a  maximum  rate,  which  is 
diaracteristic  of  the  microbial  metabolism.  Simldns  and  Alexander  (1984)  have  indicated 
that  degradation  in  a  laboratory  microcosm  can  be  expected  to  follow  Monod  kinetics  for  a 
microbial  population  not  limit^  by  nutrients  and  with  a  sufficient  substrate  concentration. 
However,  fim-oider  kinetics  are  often  observed  at  low  substrate  concentrations.  Because 
organic  solute  concentrations  in  the  experiment  were  quite  tow,  this  conclusion  is  consistent 
vHth  experimental  results.  Larson  (1979)  notes  that  the  rate  of  decrease  in  concentration  as  a 
fhnction  of  substrate  concentration  can  often  be  expressed  as  a  fust-order  equation,  and  that 
first-tnder  kinetics  are  generally  expected  for  biodegradation  at  low  organic  substrate 
concentrations. 

Towards  the  end  of  this  field  study  the  solute  degradation  rate  decreases,  and  this 
change  is  concurrent  with  reduced  substrate  concemrations  due  to  both  dilution  and 
degradation.  Rubin  et  al.  (1982)  have  demonstrated  in  laboratory  experiment  that  microbial 
degradation  rates  of  aromatic  compounds  can  be  reduced  at  low  substrate  concentrations. 
Similar  sUb^rate  limitations  are  thus  a  possible  eiqiilanation  for  the  rate  decreases  shown  in 
Figures  27  and  31. 

Table  13  shows  degradation  rate  constants  calculated  from  the  near  field  data.  Initial 
mass  and  mass  remaining  for  each  sampling  time  were  used  in  a  nonlinear  least  squares 
regression  to  fit  a  first-order  exponential  function  to  the  data  to  determine  the  first-order  rate 
constant,  k  (days  '),  for  each  compound.  The  maximum  rate  constants,  k.„„,  were  calculated 
from  the  maximum  slt^s  of  the  curves  shown  in  Figure  31 .  Rate  constants  corrected  for 
organic  solute  leakage  from  the  near  field,  k„^(c),  were  obtained  by  subtracting  the  leak  rate 
constant  k|  from  values.  The  relative  values  of  the  degradation  rates  given  in  Tables  12 
and  13  appear  reasonable,  with  o-DCB  being  tlK  most  slowly  degraded  test  chemical.  These 
findings  are  in  agreement  with  Callahan  et  al.  (1979). 
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Figure  30.  Plot  of  the  Relative  Tritium  Concentration  in  the  Near-Field  Region. 
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DAYS  AFTER  INJECTION 


Figure  31.  Plot  of  the  Organic  Compound  Masses,  Based  on  Near-Field  Data. 
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TABLE  13.  DEGRADATION  RATE  CONSTANTS  FOR  ORGANIC 
COMPOUNDS  IN  THE  NEAR  FIELD 


k  (days  ') 

(days  ') 

k^(c)  (days  ') 

benzene 

0.0079 

0.0120 

0.0066 

p-xylene 

0.0106 

0.0196 

0.0141 

naphthalene 

0.0072 

0.0118 

0.0063 

o-dichlorobenzene 

0.0060 

0.0114 

0.0059 

Rate  constants  were  based  on  the  analysis  of  stable  compounds  by  gas 
chromatography.  Ring  labeled  p-xylene  was  used  to  show  that  it  was  degraded  rather 
than  lost  in  the  aquifer.  p-xylene  was  oxidized  to  ^*C  labeled  intermediates  and  carbon 
dioxide,  which  exists  as  H2‘^C03  and  H‘^CO,'  at  the  existing  groundwater  pH  of 
approximately  5.  The  special  sampling  and  analysis  procedure  described  in  MacIntyre  et  al. 
(1993)  permitted  measurement  of  '^C  p-xylene  degra^tion  products.  '*C  in  the  carbonate 
precipitate  was  divided  by  in  the  basic  sample  (siq)ematant  plus  precipitate)  to  obtain  the 
fraction  of  p-xylene  converted  to  COj.  The  amount  of  '^C  in  the  water  after  extraction,  was 
divided  by  in  the  unextracted  water  sample  to  obtain  the  fraction  of  p-xylene  converted  to 
polar  intermediates  and  COj. 

Applicability  of  the  leaking  reactor  model  is  site  specific,  and  its  success  with  the 
MADE-2  test  data  is  a  fortuitous  circumstance,  dependent  on  the  arbitrarily  selected  positions 
of  the  injection  wells  in  the  Columbus  aquifer.  It  is  encouraging  that  the  reactor  and  whole 
field  methods  of  calculating  biodegradation  rates  gave  similar  results.  This  fact  provides 
more  confidence  in  the  descriptive  model  and  in  the  utility  of  these  rates  for  predictive 
purposes,  including  the  design  of  further  experiments  at  the  site. 

Degradation  of  the  organic  solutes  was  confirmed  by  a  special  sampling  of  wells  in 
high  and  low  solute  concentration  regions  of  the  plume  at  421  days  after  injection.  Selection 
of  the  wells  was  based  on  known  concentration  data  from  the  previous  snapshot  sampling. 
Each  water  sample  was  subsampled  and  processed  in  the  field.  A  10  mL  subsample  was 
placed  in  a  glass  ampoule,  made  basic  with  1  mL  of  3N  NaOH  and  1  mL  of  3N  Ca(N03)2 
and  immediately  flame-sealed.  A  separate  2  mL  subsample  was  mixed  with  scintillation 
cocktail  and  counted.  A  final  20  mL  subsample  was  extracted  with  2  mL  of  unlabeled 
p-xylene.  One  mL  of  the  resulting  xylene  layer  was  transferred  to  scintillation  cocktail  and 
counted  to  determine  the  amount  of  undegraded  p-xylene  which  remained  in  the  organic 
layer.  Two  mL  of  the  extracted  water  was  transferr^  to  scintillation  cocktail  and  counted. 
After  they  were  returned  to  the  laboratory,  the  ampoules  were  centrifuged  to  separate  a 
precipitate  containing  Ca.^*COj,  aiKl  then  they  were  broken  open.  The  precipitate  and  2  ml  of 
the  supernatant  water  were  placed  in  separate  scintillation  vials  with  cocktail,  and  counted. 

Conversion  results  for  p-xylene  in  well  samples  from  high  concentration  regions  of 
the  plume  are  given  in  Table  14  using  carbonate  precipitation,  and  in  Table  IS  based  on 
extraction  with  p-xylene.  In  Table  14,  DPM  is  calculated  by  correcting  the  precipitate 
counts  for  background  and  for  ‘^C  carry  over  from  the  solution.  Percent  conversion  is 
calculated  from  the  ‘^COa  DPM  /  Total  DPM  quotient. 
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TABLE  14.  p-XYLENE  CONVERSION  TO  COj  DETERMINED 
BY  THE  PRECIPITATION  OF  CARBONATE  FROM 
WATER  SAMPLES  TAKEN  AT  DAY  421 


Precipitate 

Solution  DPM/ml 

‘^COj 

Total  “C 

WeU 

Port 

’H 

'H 

DPM 

DPM 

%  Conversion 

D18 

11 

4532 

700 

2715 

28 

66 

94 

71 

D18 

13 

4827 

932 

5149 

82 

82 

164 

50 

D21 

13 

2544 

742 

2448 

16 

71 

87 

82 

D6 

12 

2479 

928 

3420 

38 

85 

123 

70 

D6 

14 

6033 

1760 

6952 

87 

164 

250 

65 

M30 

19 

3017 

841 

2036 

13 

81 

94 

87 

M3S 

11 

1389 

411 

1646 

15 

37 

52 

71 

M3S 

13 

1210 

354 

698 

10 

33 

42 

78 

M36 

13 

1081 

366 

933 

6 

34 

40 

86 

M40 

11 

515 

413 

537 

16 

38 

53 

71 

M41 

13 

3570 

1642 

3236 

46 

158 

204 

77 

M41 

15 

4697 

1145 

3363 

29 

110 

139 

79 

M41 

17 

1941 

403 

1268 

11 

37 

48 

77 

M43 

13 

5371 

1095 

4470 

41 

104 

145 

72 

M43 

15 

2079 

420 

1591 

10 

39 

49 

80 

M4S 

13 

2648 

708 

1646 

13 

68 

80 

84 

MSO 

11 

2163 

445 

1688 

14 

41 

55 

75 

M50 

13 

2201 

438 

1989 

23 

40 

63 

63 

Conversion  percentages  given  in  Table  IS  are  based  on  the  assumption  that  ’^C 
p-xylene  partitions  almost  completely  into  the  xylene  phase  and  '^COj  partitions  almost 
c<niq>letely  into  the  water  phase.  assumption  is  suf^rted  empiri^ly  by  the  agreement 
between  p>xylene  conversion  pecemages  determined  by  the  p-xylene  extraction  and  the 
carbonato  precipitation  method.  However,  the  behavior  of  the  COj  and  xylene  buffered 
water  system  cannot  be  predicted  from  existing  theory,  and  the  literature  does  not  contain 
experimental  measurements  on  this  system.  One  of  ^  authors  intends  to  make  and  publish 
these  measurements  in  the  near  future,  as  the  system  contains  some  chemical  curiosities.  For 
example,  COj  is  numy  times  more  soluble  in  pure  p-xylene  than  in  pure  water,  but  field 
measurements  here  indicated  that  most  of  the  COj  produced  by  p-xylene  degradation  went 
into  the  water  phase. 

Since  this  matter  was  deemed  worthy  of  further  investigation,  a  series  of  extraction 
tests  was  done  in  the  laboratory  to  determine  the  partititioning  of  H'^COs'  between  water  and 
p-xylene  in  a  volume  ratio  of  10:1  at  water  pH  values  of  4.8  and  7.0.  It  is  apparent  from 
the  results  given  in  Table  16,  for  five  replicate  extractions  at  each  pH,  that  most  of  the 
partititioned  into  the  water  phase,  even  at  a  pH  of  4.8.  The  ‘^C02  behavior  in  these 
extraction  tests  was  similar  to  that  measured  by  the  p-xylene  extraction  method  used  in  the 
field,  thus  confuming  the  effectiveness  of  the  field  procedure. 
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TABLE  15.  p-XYLENE  CONVERSION  TO  COj  DETERMINED  BY  XYLENE 
EXTRACTION  OF  WATER  SAMPLES  TAKEN  AT  DAY  421 


Well 

Port 

Phase 

^H/mL 

'"C/mL 

%  Conversion 

D18 

11 

Aqueous 

2978 

76 

D18 

11 

Organic 

1 

8 

84 

D18 

11 

Water 

3221 

91 

D18 

13 

Aqueous 

5682 

154 

D18 

13 

Organic 

0 

38 

77 

D18 

13 

Water 

5936 

199 

D21 

13 

Aqueous 

2640 

90 

D21 

13 

Organic 

1 

7 

94 

D21 

13 

Water 

2657 

96 

D6 

12 

Aqueous 

3583 

117 

D6 

12 

Organic 

1 

15 

86 

D6 

12 

Water 

3781 

136 

D6 

14 

Aqueous 

7793 

244 

D6 

14 

Organic 

2 

31 

87 

D6 

14 

Water 

8162 

281 

M3 

19 

Aqueous 

2304 

80 

M3 

19 

Organic 

1 

11 

84 

M3 

19 

Water 

2266 

95 

M35 

11 

Aqueous 

1611 

65 

M35 

11 

Organic 

5 

13 

85 

M35 

11 

Water 

1618 

77 

M35 

13 

Aqueous 

648 

37 

M35 

13 

Organic 

1 

6 

79 

M35 

13 

Water 

651 

46 

M36 

13 

Aqueous 

988 

41 

M36 

13 

Organic 

1 

6 

152 

M36 

13 

Water 

972 

27 
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TABLE  15.  p-XYLENE  CONVERSION  TO  CO^  DETERMINED  BY  XYLENE 
EXTRACTION  OF  WATER  SAMPLES  TAKEN  AT  DAY  421 


WsU 

Port 

Phase 

‘^C/mL 

%  Conversion 

M40 

11 

Aqueous 

296 

42 

M40 

11 

Organic 

0 

12 

71 

M40 

11 

Water 

308 

59 

M41 

13 

Aqueous 

3907 

182 

M41 

13 

Organic 

1 

40 

83 

M41 

13 

Water 

3893 

220 

M41 

15 

Aqueous 

4187 

128 

M41 

15 

Organic 

1 

15 

93 

M41 

15 

Water 

4209 

137 

M41 

17 

Aqueous 

1460 

44 

M41 

17 

Organic 

2 

5 

91 

M41 

17 

Water 

1424 

49 

M43 

13 

Aqueous 

5216 

125 

M43 

13 

Organic 

2 

27 

77 

M43 

13 

Water 

5244 

162 

M43 

15 

Aqueous 

1724 

43 

M43 

15 

Oirganic 

1 

5 

87 

M43 

15 

Water 

1720 

50 

M4S 

13 

Aqueous 

1796 

64 

M4S 

13 

Organic 

0 

11 

82 

M4S 

13 

Water 

1762 

77 

M50 

11 

Aqueous 

1961 

55 

M50 

11 

Organic 

5 

7 

87 

M50 

11 

Water 

1976 

63 

M50 

13 

Aqueous 

2218 

60 

M50 

13 

Organic 

1 

12 

77 

M50 

13 

Water 

21% 

78 

68 


TABLE  16.  LABORATORY  TEST  OF  ‘^COj  PARTITIONING 
BETWEEN  WATER  AND  p-XYLENE 


pH  =  4.8 


SAMPLE 

DPM/mL 

DPM/mL  (corrected) 

DPM  (total) 

Aqueous 

1226 

1191 

23820 

Phase 

Organic 

1904 

1869 

3737 

Phase 

Fraction  of  in  the  organic  phase  =  0.157 


pH  =  7 

SAMPLE 

DPM/mL 

DPM/mL  (corrected) 

DPM  (total) 

Aqueous 

Phase 

1607 

1572 

31446 

Organic 

Phase 

168 

133 

266 

Fraction  of  ‘^COj  in  the  organic  phase  =  0.008 


The  ‘^C  p-xylene  conversion  values  compare  well  with  mass  losses  calculated  from 
the  moments  analysis  of  the  p-xylene  remaining  in  the  plume.  Mass  loss  at  day  421  was 
obtained  by  interpolation  between  the  mass  losses  at  328  and  440  days.  It  is  assumed  that  all 
of  the  organics  injected,  because  of  their  structural  similarity,  were  degraded  and  not  lost  to 
other  processes. 

Conversion  values  in  Table  17  imply  that  most  of  the  p-xylene  went  to  energy 
production  and  was  neither  incorporated  into  new  biomass,  nor  lost  to  sorption  or 
volatilization.  Conversions  were  calculated  separately  for  high  ( >  Ippm)  and  low  ( <  Ippm) 
solute  concentration  regions  of  the  plume.  The  difference  in  degradation  between  these  two 
regions  of  the  aquifer  were  not  found  to  be  significant.  Some  of  this  difference  may  have 
been  due  to  analytical  errors  (e.g.  loss  of  Ca'^CO,  during  precipitate  transfer  prior  to 
counting,  or  COj  dissolution  in  xylene).  The  similarity  of  these  results  imply  nonlocalized, 
concentration  independent  degradation.  Toxic  effects  and  nutrient  or  oxygen  limitations  on 
the  microbiota  are  not  indicated. 
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TABLE  17.  DEGRADATION  OF  p-XYLENE  MEASURED  IN 
WATER  SAMPLES  TAKEN  AT  421  DAYS  AFTER  INJECTION, 
EXPRESSED  AS  A  WEIGHT  PERCENT 

Converted  to  CO2  Converted  to  CO2  Mass  Loss  From 
(precipitation  metluxl)  (  extraction  method)  moments  analysis 

[X]>lppm(n=8)  72.7  (±11.5)  85.1  (±  6.3)  98 

[Xl<lppm(n=10)  75.7  (±6.9)  82.6  (±6.1)  98 


These  radiochemical  results  for  '^C  p-xylene  degradation  to  '^C02  provided 
quantitative  support  for  a  qualitative  observation  made  in  Section  5  regarding  the  spatial 
distributions  of  p-xylene  and  of  ‘^C  in  the  snapshot  samples.  The  results  indicated  that  the 
bulk  of  the  p-xylene  remained  in  the  near  rield  throughout  the  study,  while  ‘^C  concentrations 
increased  signiricantly  in  the  far  field,  and  even  appeared  to  be  transported  along  with  the 
tritium  in  the  far  rield.  These  distributions  can  o^y  be  explained  by  the  degradation  of 
p-xylene.  The  ^*C  p-xylene  had  to  have  been  converted  to  ‘^COj  in  the  near  field  through 
chonical  and/or  biological  processes,  and  was  then  transported  to  the  far  field  with  the 
tritium  tracer.  The  stability,  or  conservative  nature,  of  '*C02  in  the  Columbus  aquifer  was 
mainly  due  to  the  acidic  groundwater  conditions,  which  prevented  the  formation  of  carbonate 
mineials.  Equilibrium  thermodynamic  model  calculations  showed  that  the  grou  idwater  is 
highly  undersaturated  with  respect  to  calcite.  Ihe  absence  of  carbonate  minerals  inferred 
from  these  calculations  was  confirmed  by  the  similarity  of  total  carbon  and  organic  carbon 
contents  (at  the  0.(X)5  percent  level)  determined  using  a  Leco  Corp.  (St  Joseph,  MI)  Model 
WR112  carbon  analyzer.  Investigators  considering  use  of  the  above  radiochemical  method 
should  first  determine  whether  the  chemistry  of  the  aquifer  of  concern  will  permit  its  use 
effectively.  Conditions  which  would  promote  the  precipitation  of  carbonate  minerals  in  the 
saturated  zone  of  an  aquifer  would  preclude  the  use  of  the  above  method  for  measuring 
the  degradation  of  organic  compouiids. 

TIME-SERIES  DATA  ANALYSIS 

Concentration  time-series  data  for  several  sampling  points  in  the  near  field  were  used 
to  estimate  the  sorption  of  organic  compounds  during  the  field  study,  and  to  provide 
ind^)endent  computations  of  the  organic  degradation  rates  for  comparison  with  those  derived 
ficm  the  qntial  moments  analysis.  Retardation  factors  and  degradation  rates  were  estimated 
Ity  fitting  the  time-series  data  to  an  analytical  solution  of  the  one-dimensional  convection- 
dispersmn  equation  coittaining  linear  equilibrium  sorption  and  first-order  decay  terms. 

Normalized  solute  breakthrough  curves  (BTCs)  for  three  sampling  points  in  the  near 
field  are  shown  on  Figure  32.  These  data  were  typical  for  the  site  in  the  following  respects. 
First,  none  of  the  BTCs  showed  complete  passage  of  the  tritium  pulse.  The  slow  transport 
of  the  tracer  fimn  the  low  hydraulic  conductivity  zone  in  the  near  field  precluded  complete 
ratmittmiqt  of  tracer  breakthrough.  Second,  the  normalized  tritium  concentrations  for  each 
sampUng  point  were  consistently  the  highest  of  all  the  solutes,  as  was  expected.  However,  a 


Figure  32,  Solute  Breakthrough  Curves  at  Selected  Points 
(x.  y.  z  coordinates  are  in  meters) 
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yijRial  UMjiectkm  of  the  breakthrough  responses  indicates  that  the  area  under  the  tritium  BTC 
tyas  cootisteiitly  larger  duui  those  of  all  of  the  other  solutes.  One  would  expect,  based  on 
chronuittqpii^rtuc  theoiy,  diat  the  areas  under  all  tracer  BTCs  would  be  ai^roximately  e^ial, 
a«iiiipi;iy  (1)  omqjlete  or  nearly-complete  breakthrough  of  all  solutes,  (2)  the  absence  of 
inoceaa^  wbich  remove  tracer  from  solution  (e.g..  irreversible  adsorption  or 
faiotr^U^t^naation),  and  (3)  steady-flow  conditions.  For  the  organic  con^)ounds  and 
areas  un^  die  normalized  BTCs  are  smaller  for  benzene,  naphthalene,  and  p-xylene  than  for 
o-DCB  and  These  results  are  generally  consistent  with  the  ^latial  moments  results, 
which  indiotted  die  lowest  mass  recoveries  from  the  most  biodegradable  organics  (i.e., 
boszene,  nqdithalene,  and  p-xylene)  with  consistendy  higher  recoveries  for  o-DCB. 

PEC^LADATION  RATES  AND  RETARDATION  FACTORS 


Due  to  the  geohydrologic  conditions  in  the  Columbus  aquifer,  retardartion  facttns  for 
the  organk;  solutes  cannot  be  reliably  estimated  from  a  conventional  comparison  of 
corresponding  first  tenqioral  moments  for  the  organic  solutes  and  for  tritium.  Inflection  of 
^  poiig  tinMi-'Seiies  data  indicates  that  the  mean  arrival  times  of  the  organic  confounds  are 
less  Aan  ^  correfwnding  mean  arrival  times  of  tritium  in  most  cases,  leading  to  obviously 
emxcous  retardation  factors  of  less  than  one.  The  source  of  this  apparent  error  is  that 
convendonal  mean  arrival  time  analysis  does  not  account  for  the  degradation  of  organic 
compounds,  as  foffod  during  this  field  study,  and  thus  cannot  yield  meaningful  retardation 
factors.  An  alterp<ative  procedure  was  used  to  calculate  retard^ion  factors  for  the  organic 
solutes  that  incluijled  solute  degradation.  In  this  procedure,  the  time-series  concentration  data 
for  tritipm  and  organic  solutes  were  fitted  to  an  analytical  solution  of  the  one-dimensional 
convectitm-difiersion  equation  which  included  both  linear  equilibrium  sorption  and  first-order 
decay  terms. 

One-dimensional  steady  tranfiort  of  a  pulse  release  of  solute  undergoing  linear 
equilitenun  adsorption  and  first-order  decay  is  described  by: 


-  n^'  -  1/:^  ^  _  ac 

dt  ^  dx  ti  dt  ^  ' 


(21) 


In  the  case  modeling  this  field  experiment,  the  source  boundary  condition  is. 


Cb,  0  <  f  s  0 

^0.  f>^,  J 


(22) 


where  c^  is  die  vobime-avoaged  resident  conceittration  of  the  solute  in  the  liquid  phase,  s  is 
die  sorted  concentration  per  unit  mass  of  the  solkl  iduse,  x  is  distance,  D  is  a  dispersion 
coefficient  reflecting  die  condiined  effects  of  diffiiskm  and  hydrodynamic  difierskm,  V  is  the 
average  p(»e-water  velocity,  is  the  bulk  mass  density,  i;  is  the  porosity,  and  p  is  die  first- 
order  detey  rale  craostant. 

A  nofdinear  least-squares  inversion  code  develt^ied  by  Parker  and  van  Genuchten 
(|98#)  was  used  to  estimate  the  solute  transport  paramraers  associated  with  Equations  (21) 
and  (22),  fiom  the  time-series  concentration  data  for  discrete  sampling  points  located  in  the 
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near  field.  In  general,  the  procedure  consisted  of  fitting  the  tritium  breakthrough  response 
for  a  specific  sample  point  to  obtain  estimates  for  the  physical  transport  parameters  V  and  D. 
In  estimating  R  and  n  for  each  organic  compound  from  their  breakthrough  curves,  these  V 
and  D  values  v^ere  used  as  constants,  and  a  fitting  procedure  was  applied  to  obtain  values  for 
R  and  fi.  Time-series  data  for  11  BarCad*  samplers  and  6  MLS  points  were  selected  for 
analysis.  Only  data  for  sampling  points  exhibiting  complete  (or  nearly  complete)  single-peak 
breakthrough  responses  were  analyzed.  Breakthrough  responses  exhibiting  multiple  peaks 
were  common  among  the  sampling  points  evaluated,  but  were  excluded  from  the  analysis 
because  they  did  not  satisfy  tte  assumptions  inherent  in  Equation  21. 

Physical  transport  parameter  estimates  for  the  17  discrete  sampling  points  are 
presented  in  Table  18.  In  addition,  statistical  summaries  of  the  retar^tion  factors  and  rate 
constants  for  the  organic  compounds  are  given  in  Tables  19  and  20.  The  mean  retardation 
factors  for  benzene,  naphthalene,  p-xylene,  and  o-DCB  are  1.20,  1.45,  1.16,  and  1.33, 
respectively.  The  mean  first-order  rate  constants  describing  degradation  of  benzene, 
naphthalene,  p-xylene,  and  o-DCB  are  0.010,  0.013,  0.016,  and  O.OOS  d  ',  respectively. 

Note  that  the  estimates  of  D  presented  in  Table  18  for  the  tritium  data  may  not  be 
reliable.  Several  features  of  this  natural  gradient  experiment  (e.g.,  extreme  aquifer 
heterogeneiQr,  three-dimensional  non-steady  flow,  a^  a  finite-volume  tracer  solution  source) 
violate  the  simplistic  physical  transport  assumptions  embodied  in  Equation  (21).  The 
reported  estimates  of  V  and  D  represent  a  set  of  effective  parameters  which  together  produce 
a  physical  tranport  response  that  approximates  the  tritium  breakthrough  data.  Our  primary 
interest  here  is  in  the  difference  between  the  tritium  and  organic  solute  breakthrough 
responses  caused  by  sorption  and  degradation,  and  not  estimation  of  dispersion  coefficients. 

Transformation  rates  for  the  organic  compounds  were  independently  estimated  from 
the  solute  mass  balances  (with  results  given  in  Table  12)  and  from  the  time-series  data.  A 
first-order  rate  process  was  assumed  to  approximate  the  overall  degradation  process  in  both 
analyses.  A  comparison  of  the  estimated  rate  constants  obtained  by  the  two  methods  is  given 
in  Table  21,  and  indicates  little  difference  between  the  methods  for  benzene,  p-xylene,  and 
o-DCB. 

The  disaipearance  and  transformation  of  the  organic  solutes  during  this  experiment 
demonstrated  that  natural  degradation  processes  were  able  to  effectively  reduce  th^  levels 
of  dissolved  organic  contaminants  in  a  reasonable  time  frame.  This  result  suggests  that,  for 
similar  organic  conpounds,  it  might  be  best  to  restrict  aquifer  remediation  activities  to  the 
contaminant  source  region.  Possible  techniqiKs  could  include  pumping  to  remove  NAPLs 
(non-aqueous  (diase  liquids),  or  excavation  of  contaminated  aquifer  material.  The  source 
would  then  be  reduced,  but  probably  would  twt  be  elimiiuited  entirely  by  this  treatment.  In 
an  aquifer  with  approximately  steady  flow,  and  the  correct  physical,  chemical  and  biological 
conditions,  the  plume  of  organic  solutes  from  the  reduced  source  could  attain  a  steady-state 
limit.  In  such  an  apiifer,  the  spatial  boundary  of  this  plume  would  be  determined  primarily 
by  the  hydrology  and  geochemistry  of  the  site,  solute  sorption,  biodegradation  by  indigenous 
mictobes.  redox  opacity,  arid  oxygen  and  nutrient  supply. 

This  observation  represents  the  most  important  result  of  the  field  study,  and  could 
have  inportant  inplications  for  the  restoration  of  aquifers  contaminated  by  organic 
chemicals.  Active  remediation  would  not  be  needed  in  situations  where  natural  degradation 
rates  were  sufficient  to  reduce  contaminant  concentrations  to  safe  levels  before  they  were 
tranported  to  regions  presenting  risks  to  human  health  or  the  environnwnt.  Remediation  by 
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TABLE  1«.  VELOCITY  AND  DISPERSION  COEFFICIENT  ESTIMATES 
FROM  TRITIUM  TIME-SERIES  DATA 


INsperdoo 

Stiwnlir 

Velocity, V 

CoeffIdeiityD 

Point 

x(m) 

y(in) 

*(m) 

(m/d) 

(mVd) 

r* 

B61B02 

-4.0 

8.7 

60.03 

0.0433 

0.0075 

0.80 

B62A02 

-0.7 

7.2 

59.71 

0.0227 

0.0230 

0.98 

B62B01 

-2.4 

7.9 

56.66 

0.0687 

0.2303 

0.87 

B63AQ2 

2.6 

6.1 

57.02 

0.0412 

0.1728 

0.95 

B63B01 

1.1 

6.8 

55.46 

0.1531 

0.8392 

0.96 

B64A01 

5.5 

5.0 

55.08 

0.1098 

0.3497 

0.96 

B64AQ2 

5.5 

5.0 

59.66 

0.0601 

0.0497 

0.90 

B68A02 

-2.5 

18.0 

58.12 

0.1187 

0.2438 

0.90 

B68BQ2 

0.1 

17.0 

59.98 

0.0859 

0.0769 

0.79 

B69A02 

2.1 

16.3 

58.53 

0.0533 

0.0823 

0.98 

B70BQ2 

5.3 

15.4 

60.01 

0.0668 

0.0330 

0.94 

M06S1S 

-0.5 

19.1 

59.85 

0.0591 

0.5109 

0.87 

M06Sn 

-0.5 

19.1 

60.61 

0.1020 

0.1551 

0.94 

M06615 

3.4 

17.7 

59.53 

0.0094 

0.2893 

0.98 

M06617 

3.4 

17.7 

60.29 

0.1221 

0.2400 

094 

M066i9 

3.4 

17.7 

61.05 

0.1039 

0.1614 

0  98 

MQ6719 

7.1 

16.4 

61.04 

0.1265 

0.4454 

0.92 

n 

17 

17 

mean 

0.0792 

0.2300 

std.  dev. 

0.0389 

0.2082 

95%  Cl 

0.0200 

0.1071 
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TABLE  19.  RETARDATION  FACTOR  ESTIMATES  FROM 
HYDROCARBON  TIME-SERIES  DATA 


Sample 

Retardation  Factor.  R 

Pimt 

Benzene 

Naphthalene 

p-Xylene 

o-DCB 

B61B02 

1.00 

1.00 

1.02 

1.08 

B62A02 

1.08 

1.26 

1.02 

1.14 

B62B01 

1.15 

1.61 

1.20 

1.14 

B63A02 

1.00 

1.56 

1.00 

1.40 

B63B01 

1.71 

2.00 

1.46 

1.54 

B64A01 

— 

— 

— 

1.41 

B64A02 

1.00 

1.11 

1.00 

1.09 

B68A02 

1.06 

1.21 

1.08 

1.06 

B68B02 

1.24 

1.41 

1.14 

1.18 

B69A02 

1.29 

1.57 

1.02 

1.36 

B70B02 

~ 

— 

— 

1.23 

M06515 

— 

— 

1.47 

M06517 

— 

1.12 

M0661S 

1.32 

2.00 

1.28 

1.77 

M06617 

1.14 

1.23 

1.13 

1.08 

M06619 

1.36 

— 

1.69 

1.65 

M06719 

— 

1.00 

1.84 

n 

12 

11 

13 

17 

mean 

1.20 

1.45 

1.16 

1.33 

St.  dev. 

0.20 

0.32 

0.20 

0.25 

95%  Cl 

0.13 

0.21 

0.12 

0.13 
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TABLE  20.  ESTIMATES  OF  FIRST-ORDER  BIOTRANSFORMATION 
RATE  CONSTANTS  FROM  TIME-SERIES  DATA 


goHUpiy 

Point 

Bouene 

Naphthaioie 

p-Xylrae 

o-DCB 

B61B02 

0.0113 

0.0120 

0.0148 

0.0035 

B62AQ2 

0.0024 

0.0022 

0.0029 

0.0005 

B62B01 

0.0160 

0.0214 

0.0138 

0.0046 

B63A02 

0.0144 

0.0257 

0.0211 

0.0132 

B63B01 

0.0104 

0.0163 

0.0053 

0.0027 

B64A01 

— 

— 

— 

0.0118 

B64A02 

0.0150 

0.0183 

0.0151 

0.0055 

B68A02 

0.0132 

0.0107 

0.0094 

0.0027 

B68BQ2 

0.0051 

0.0090 

0.0118 

0.0038 

B69A02 

0.0070 

0.0083 

0.0124 

0.0045 

B70B02 

— 

— 

— 

0.0049 

M0651S 

— 

— 

— 

0.0037 

M065n 

— 

— 

— 

0.0056 

M0661S 

0.0039 

0.0058 

0.0155 

0.0022 

M06617 

0.0094 

0.0095 

0.0091 

0.0028 

M06619 

0.0105 

— 

0.0164 

0.0103 

M06719 

— 

— 

0.0547 

0.0094 

n 

12 

11 

13 

17 

mean 

0.0099 

0.0127 

0.0156 

0.0054 

St.  tlev. 

0.0043 

0.0067 

0.0122 

0.0035 

95%  a 

0.0027 

0.0045 

0.0074 

0.0018 
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natural  degradation  processes  is  likely  to  be  applicable  to  residual  organic  contaminants  (i.e., 
NAPL  organics  immobilized  by  capillary  forces  that  cannot  be  effectively  removed  from  the 
aquifer  by  pumping)-  The  monetary  and  environmental  cost  savings  of  allowing  natural 
biological  restoration  of  the  residual  contaminants  are  potentially  enormous. 

The  degradation  rates  determined  for  the  four  aromatic  compounds  in  the  Columbus 
aquifer  are  now  being  used  in  the  design  and  modeling  stages  of  a  new  field  test  which  will 
use  an  enqilaced  NAPL  source.  This  test  at  the  Columbus  site  will  be  a  demonstration  of  the 
practicality  of  natural  degradation  (natural  attenuation)  as  a  remediation  action  for  a  steadily 
leaching  source  in  an  aquifer.  Assuming  that  a  steady-state  contaminant  concentration 
situation  is  attained  in  this  test,  natural  degradation  will  be  a  verified  groundwater 
contaminant  treatment  option.  It  is  proposed  that  experiments  similar  to  this  one  should  be 
done  at  several  other  test  sites  in  order  to  prove  the  general  applicability  of  this  method,  and 
for  determining  accurate  aerobic  degradation  rates  for  use  in  the  development  of  more 
accurate  fate  and  transport  models.  The  test  sites  could  then  be  used  to  conducted  NAPL 
tests  which  might  further  support  the  natural  attenuation  treatment  option,  and  might  even 
suggest  this  technique  as  the  action  of  Hrst  choice  in  certain  groundwater  contamination 
situations. 


TABLE  21.  SUMMARY  OF  FIRST-ORDER  DEGRADATION  RATE 
CONSTANTS  AND  HALF-LIVES  OF  AROMATIC  COMPOUNDS 


Tracer 

Method 

Rate  Constant  (d  ') 

Half-Life  (d) 

benzene 

SMA 

0.007 

87 

TA 

0.010 

69 

naphthalene 

SMA 

0.006 

99 

TA 

0.013 

53 

p-xylene 

SMA 

0.011 

63 

TA 

0.016 

43 

o-DCB 

SMA 

0.005 

116 

TA 

0.005 

139 

SMA  =  spatial  moments  analysis 
TA  =  time-series  analysis 


The  rate  constants  estimated  by  the  two  methods  for  naphthalene  are  significantly 
different,  i.e.,  rate  estimates  differ  by  a  factor  of  about  two.  The  explanation  for  the 
difference  between  the  naphthalene  rate  constants  is  not  evident.  In  general,  results  for  both 
methods  are  consistent  with  the  expected  relative  biodegradabilities  of  the  four  compounds. 
The  lowest  transformation  rates  are  estimated  for  o-DCB,  the  least  degradable  organic,  while 
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l^glier  rates  are  indicated  for  benzene,  p-xylene,  and  ni4}hthaiene,  which  have  been  slwwn  in 
Mwnttory  studies  to  be  more  susceptible  to  biodegradation  (Tabak  et  al.,  1981). 

COMPARISON  OF  SORPTION  IN  LABORATORY  AND  FIELD  EXPERIMENTS 

As  shown  in  Table  22.  the  mean  retardation  factors  estimated  for  these  organic 
con^punds  from  the  time-series  data  for  discrete  sampling  points,  are  remarkably  consign! 
with  die  laboratory  measured  values  reported  by  MacIntyre  et  al.  (1993)  using  aquifer 
nuderial  sanqiles  from  the  Columbus  site.  Furthermore,  the  same  relationship  in  the  relative 
magnitudes  of  the  retardation  factors  among  the  four  organic  compounds  is  indicated  for  both 
the  field  and  labortfory  results. 


TABLE  22.  COMPARISON  OF  AVERAGE  RETARDATION  FACTORS 
ESTIMATED  FROM  LABORATORY  TESTS  AND  FIELD  DATA 


Lab 

Field 

benzene 

1.30 

1.20 

naphthalene 

1.42 

1.45 

p-xylene 

1.24 

1.16 

o-DCB 

1.32 

1.33 

The  mecluninn  reqxmsible  for  the  sorption  of  these  compounds  is  unknown.  The 
conqxwnds  range  from  sli^tly  to  moderately  hydrc^hobic,  as  indicated  by  their  octanol- 
water  partitioning  coefficients:  benzene»13S;  naphthalene =20(X);  p-xylene=1410: 
o-DCB»2S()0  (Oherini  et  al.,  1989).  For  soils  ^  aquifer  materials  with  organic  carbon 
CQittents  of  greater  than  0. 1  percent,  the  sorption  mechanism  for  hydrophobic  organics  is 
assumed  to  involve  the  partitionii^  of  the  sorbate  into  organic  matter  tMt  is  bound  to  the 
mineral  surfaces  (Kanckhoff  et  al.,  1979).  However,  the  average  organic  carbon  content  of 
this  aquifer  is  tnuch  less  than  0.1  percent  (MacIntyre  et  al.,  1991).  As  might  be  expected 
from  such  a  low  carbon  content  material,  Macintjm  et  al.  (1991)  found  a  low  correlation 
(r^aO.048)  between  the  K.,  for  naphthalene  and  organic  carbon  content  in  their  analysis  of  SO 
a^piifer  samples  collected  from  the  test  site.  Although  these  investigators  did  not  identify  the 
naphthalene  sorption  mechanism,  they  did  conflrm  through  a  series  of  batch  sorption  kinetic 
and  equilibrium  eiqieriments  that  (1)  the  sorption  process  was  rapid  with  equilibrium 
achieved  within  24  hours,  and  (2)  the  naphthalene  isotheim  was  linear  in  tte  range  of 
OrS  mg/L.  These  findings  sui^rt  our  presumption  of  linear  equilibrium  sorption  reactions 
in  the  parameter  estimation  aiialysis  of  the  time-series  data  for  the  organic  compminds. 

The  similarity  between  the  field  and  laboratory  retardation  factors  is  enctniraging.  It 
suggests  that  anall-scale  laboratory  tests  can  provide  meaningful  estimates  of  field  average 
sorption  inx^ierties  when  measurements  are  sufficient  to  capture  the  Held  variability. 
However,  we  qualify  this  by  noting  that  transport  scale  of  the  time-series  data  was  only  on 
die  order  of  20  meters.  This  scale  of  observation  is  considerably  smaller  than  the  transport 
scales  usually  connected  with  groundwater  contamination  problems  (i.e.,  tens  to  hundred  of 
meters). 
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SECTION  Vll 
AQUIFER  DISPERSIVITY 


Dispersivity  in  the  aquifer  can  be  calculated  from  the  first  and  second  spatial  moments 
of  the  tritium  plume  presented  in  Section  5.  Because  the  mean  velocity  field  at  the 
Columbus  site  is  nonuniform,  dispersivity  cannot  be  predicted  directly  from  the  rate  of 
change  of  the  second  moments  as  has  been  done  for  other  experiments  at  other,  more 
homogeneous  sites  (e.g.,  Freyberg,  1986;  Garabedian  et  al.,  1991)  Adams  and  Gelhar 
(1992)  developed  a  method  to  calculate  dispersivity  at  the  Columbus  site,  and  applied  it  to 
the  results  of  a  bromide  tracer  field  test.  They  derived  general  moments  equations  for  two- 
dimensional  advective-dispersive  transport  in  a  spatially  nonuniform  flow  Held,  following  an 
instantaneous  discharge  of  mass.  The  two-dimensional  velocity  field  is  defined  by. 


Vo  (  1  -I-  ay  ) 

(23) 

-Vo  ax 

(24) 

where  Vi.  is  the  longitudinal  velocity,  a  is  a  constant,  is  the  initial  longimdinal  velocity, 

Vj  is  the  transverse  velocity,  and  y  and  x  are  spatial  coordinates  aligned  with  the  longitudinal 
and  transverse  flow  directions,  respectively.  Substitution  of  this  function  into  the  general 
moments  equation  (Equation  S-S  of  Adams  and  Gelhar,  1992),  and  solution  of  the  resulting 
ordinary  differential  equations  corresponding  to  the  first  and  second  moments  (Equations  5-8, 
5-9,  and  5-12  of  Adanis  and  Gelhar,  1992)  yields  analytical  expressions  for  the  mean 
displacement,  y 

y\i)  =  c,  (exp(C2t)  -  1)  (25) 

the  longitudinal  second  moment, 


®ii(0  “  + 


e 


2(V 


Cz 


the  horiBHhtal  transverse  second  moment, 


(26) 


2A22C,  A22(Vo-C,C5j)) 


0 


(27) 


where  c,=(l  +  oiAu)/a,  C2=aVo,  A„  is  the  longitudinal  dispersivity,  and  A22  is  the 
transverse  dispersivity. 

Results  of  this  dispersion  analysis  are  summarized  in  the  four  graphs  shown  on  Figure 
33.  The  predicted  and  observed  tritium  plume  mean  displacements  over  time  are  shown  on 
Figure  33(a).  The  corresponding  longitudinal  and  horizontal  transverse  second  moments  are 
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Figure  33.  Comparison  of  (a)  Mean  Displacement,  (b)  Longitudinal  Second  Moment,  and 

(c)  Transverse  Second  Moment  From  Tritium  Data  and  Nonuniform  Flow  Model 
Open  and  Solid  Circles  Represent  Tritium  Data  With  and  Without  Extrapolation 
at  328  Days,  Respectively.  Open  Square  is  Bromide  Data  From  Adams  and 
Gelhar  (1992).  Dotted  and  Solid  Lines  Denote  Model  Predicted  Fits  to  Tritium 
Data  With  and  Without  Extrapolation.  Diagram  (d)  Shows  Mean  Longitudinal 
Velocity  as  a  Function  of  Displacement  for  Base  and  Extrapolated  Cases. 
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plotted  against  mean  displacement  on  Figures  33(b)  and  33(c).  The  bromide  plume 
displacement  and  moment  results  of  Adams  and  Gelhar  (1992)  are  also  presented  in  these 
figures  for  comparison  with  the  tritium  plume  behavior  seen  in  this  experiment.  Figure 
33(d)  shows  the  velocity  in  the  longitudinal  direction  as  a  function  of  displacement.  Two 
different  tritium  moments  are  presented  for  the  snapshot  at  328  days;  one  for  the  moments 
estimated  directly  from  the  tritium  measurements,  and  the  other  for  moments  estimated  by 
extrapolation  of  the  tritium  plume  front  to  juxount  for  it’s  truncation.  The  extrapolation 
assui^  that  the  tritium  activity  deficit  of  approximately  23  percent  observed  at  328  days 
was  entirely  due  to  truncation.  The  missing  tritium  activity  was  assumed  to  be  linearly 
distributed  within  a  region  equivalent  to  the  advance  of  the  tritium  plume  between  132  and 
224  days. 

The  parameters  assumed  or  estimated  in  fitting  Equations  (25),  (26)  and  (27)  to  the 
tritium  plume  data  are  given  in  Table  23.  Note  that,  for  both  the  base  case  and  the 
extrapolation  case,  was  assumed  to  be  0.014  m/d  (S  m/yr)  based  on  the  estimated  rate  of 
movement  of  the  plume  in  the  near  field  (Section  5).  Figure  33(a)  indicates  that  the 
longitudinal  first  moments  of  the  plume  data  are  approximated  reasonably  well  by  Equation 
(25)  for  both  cases.  The  linear  velocity  trends  in  the  longitudinal  direction  are  shown  in 
Figure  33(b).  The  longimdinal  velocity  at  a  displacement  of  1(X)  meters  is  0.46  m/d  for  the 
base  case,  and  0.79  m/d  for  the  extrapolation  case.  Both  estimates  are  below  the  seepage 
velocity  of  ^1.1  m/d  estimated  from  the  advance  of  the  tritium  front  between  132  ai^  224 
days,  and  are  therefore  within  the  realm  of  possible  mean  velocities. 


TABLE  23.  il.OW  AND  DlSPERSrVlTY  PARAMETER  ESTIMATES 

FOR  THE  TRITIUM  PLUME 


Parameter 

Base  Case 

Extrapolation  Case 

U„  (m/d) 

0.0137 

0.0137 

a 

0.328 

0.569 

A,,  (m) 

19.6 

9.5 

A22  (m) 

2.2 

— 

ff„^(0)  (m^) 

7.0 

7.0 

<T22^(0)  (m^) 

7.0 

7.0 

The  pr^icted  trends  in  the  longitudinal  second  moments  are  nearly  coincident  for  the 
base  case  and  the  extrapolation  case,  and  generally  capture  the  overall  trend  indicated  by  the 
tritium  data.  Longitudinal  dispersivities  of  9.5  and  19.6  meters  are  estimated  for  the 
extr^lation  and  base  cases,  respectively.  The  insensitivity  of  the  longitudinal  second 
moments  to  the  magnitude  of  dispersivity  over  the  expected  range  of  and  a  reflects  the 
dominant  influeixte  of  advective  transport  in  the  evolution  of  the  tritium  plume.  Because 
frontal  truncation  of  the  tritium  plume  was  evident  in  the  data,  the  extrapolation  case 
longitudinal  dispersivity  estimate  of  approximately  10  meters  is  consider^  the  more  reliable 
of  the  two  estimates. 
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A  liorisQiitMl  tnuisverse  dispersivity  of  2.2  meters  produced  the  best  agreement  with 
IQie  obierved  transverse  second  moments.  However,  this  value  is  probably  overestimated 
bnPiMe  hin^  model  used  in  the  analysis  assumes  a  convergent  flow  fieid,  when  in 
jl^  r^km  of  convogence  is  limited  to  the  southern  one-third  to  one-half  of  the  sanq)Iing 
(Ine  Fipm  3).  Cmuequently,  an  unrealistically  large  transverse  di^rsiviQr  is 
retiuiiied  in  luting  Equation  (27)  to  the  observed  transverse  second  moments,  to  compensate 
for  die  ocmvogiiig  flow  assumed  to  exist  over  the  full  length  of  the  flow  field. 

The  relative  quality  of  dispersivity  values  obtained  from  the  bromide  tracer 
exiwrinient  (Adams  and  (Selhar,  1992)  and  from  this  tritium  tracer  experiment  is  not  treated 
Ime.  ^owevo*,  trkiated  water  is  the  better  tracer,  since  solute  chemical  properties  cannot 
affect  qdcu^ated  water  transmit  prq)erties.  Also,  the  tritiated  water  smqishots  included  more 
of  field  site  was  the  case  in  sn^hots  of  the  earlier  bromide  experiment. 
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i  SECTION  vm 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  dis2y^)earaiice  and  transformation  of  the  organic  solutes  during  this  experiment 
demonstrated  that  natural  degradation  processes  were  able  to  effectively  reduce  these  levels 
of  dissolved  organic  contaminants  in  a  reasonable  time  frame.  This  result  suggests  that,  for 
similar  organic  compounds,  it  might  be  best  to  restrict  aquifer  remediation  activities  to  the 
contaminant  source  region.  By  reducing  this  source,  the  resulting  plume  of  organic  solutes 
could  possibly  maintain  a  steady-state  limit,  given  the  correct  physical,  chemical  and 
biological  conditions.  In  such  an  aquifer,  the  spatial  boundary  of  this  plume  would  be 
determined  primarily  by  the  hydrology  and  geochemistry  of  the  site,  solute  sorption, 
biodegradation  by  i^igenous  microbes,  redox  capacity,  and  oxygen  and  nutrient  supply. 

This  observation  represents  the  most  important  result  of  this  field  study,  and  could 
have  important  implications  for  the  restoration  of  aquifers  contaminated  by  organic 
chemicals.  Active  remediation  would  not  be  needed  in  situations  where  natural  degradation 
rates  were  sufficient  to  reduce  contaminant  concentrations  to  safe  levels  before  they  were 
transported  to  regions  where  they  might  be  dangerous  to  human  health  or  the  environment. 
Remediation  by  natural  degradation  processes  is  likely  to  be  applicable  to  residual  organic 
contaminants  (i.e.,  NAPL  organics  immobilized  by  capillary  forces  that  cannot  be  effectively 
removed  from  the  aquifer  by  pumping).  The  monetary  and  environmental  cost  savings  of 
allowing  natural  biological  restoration  of  the  residual  contaminants  are  potentially  enormous. 
After  a  review  of  the  current  scientific  literature,  it  appears  that  this  study  is  the  first  field 
experiment  to  prove  conclusively  that  hydrocarbon  solute  losses  were  due  to  chemical 
(tegradation  rather  than  physical  losses,  sutdi  as  sorption  and  vaporization. 

Further  research  is  needed  to  build  upon  results  of  this  study,  to  further  siq>port  the 
observation  of  natural  attenuation  and  to  attempt  to  quantify  the  environn^ntal  requirements 
for  this  process.  The  degradation  rates  determined  for  the  four  aromatic  compounds  in  the 
Columbus  aquifer  are  now  being  used  in  the  design  and  modelii^  stages  of  a  new  field  test 
which  will  use  an  enqilaced  NAPL  source.  This  test  at  the  Columbus  site  will  demonstrate 
the  practicality  of  natural  degradation  (natural  attenuation)  as  a  remediation  action  for  a 
steadily  leaching  source  in  an  aquifer.  This  experiment  more  closely  emulates  the  effects  of 
a  fuel  spill.  Assuming  that  a  steady-state  contaminant  concentration  situation  is  attained  in 
this  test,  natural  degradation  will  be  a  verified  groundwater  contaminant  treatmeitt  option. 

Experiments  similar  to  this  one  should  be  done  at  several  other  test  sites  to  prove  the 
general  ai^licability  of  this  method,  and  to  measure  both  aerobic  and  anaerobic  degradation 
rates  for  use  in  the  development  of  more  «xurate  fate  and  transport  models.  The  technique 
of  using  '^-labeled  compounds  is  suggested  for  future  field  experiments,  because  it 
distinguishes  solute  degradation  from  solute  losses  by  sorption  and  evaporation,  and  permits 
the  demonstration  of  more  accurate  mass  balances  throughout  the  course  of  the  investigation. 
If  multiple  test  sites  become  available,  they  could  then  be  used  to  conduct  NAPL  studies 
whkh  might  further  support  the  natural  attenuation  treatment  option,  and  might  even  suggest 
this  technique  as  the  action  of  first  choice  in  (Attain  groundwater  contamination  situations. 
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SECTION  X 
APPENDIX  A 


DISSOLVED  OXYGEN  MEASUREMENTS 

Dissolved  oxygen  (DO)  data  for  six  sampling  dates  are  given  in  Tables  A-1  through  A-6. 
Each  table  provides  the  MLS  or  piezometer  number,  the  measured  dissolved  oxygen  concentration 
in  mg/L,  and  the  ten^)erature  in  degrees  centigrade  for  each  water  sample  at  the  time  DO  was 
OMasured.  The  spatial  coordinates  of  each  DO  monitoring  well  are  presented  in  Table  A-7. 
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TAH4I  A-1.  LOCATIONS  OF  DISSOLVED  OXYGEN  MONITORING  SITES 


mM 

X 

(■0 

Y 

(m) 

Z 

(m) 

MLS 

Port 

X 

(m) 

Y 

(m) 

Z 

(m) 

D006 

6 

-1.5 

6.6 

56.62 

M172 

1 

3.1 

98.6 

53.85 

0006 

10 

-1.5 

6.6 

57.63 

Ml  72 

10 

3.1 

98.6 

57.28 

D006 

20 

-1.5 

6.6 

60.17 

M172 

20 

3.1 

98.6 

61.09 

0010 

3 

-2.8 

-2.3 

55.17 

M181 

1 

-1.6 

117.0 

53.80 

DOlO 

10 

-2.8 

-2.3 

57.83 

M181 

10 

-1.6 

117.0 

57.23 

DOlO 

20 

-2.8 

-2.3 

61.64 

M181 

20 

-1.6 

117.0 

61.04 

Ii04l 

1 

-.2 

3.6 

54.27 

M190 

1 

3.1 

132.1 

53.80 

m\ 

MQ41 

10 

-.2 

3.6 

57.70 

M190 

10 

3.1 

132.1 

57.23 

20 

-.2 

3.6 

61.51 

M190 

20 

3.1 

132.1 

61.04 

M066 

1 

3.4 

17.7 

54.19 

M202 

1 

1.4 

166.6 

53.65 

10 

3.4 

17.7 

57.62 

M202 

10 

1.4 

166.6 

57.07 

20 

3.4 

17.7 

61.43 

M202 

20 

1.4 

166.6 

60.88 

M07^ 

1 

2.9 

23.0 

54.07 

M235 

1 

4.4 

199.1 

53.66 

MQ7S 

10 

2.9 

23.0 

57.50 

M235 

10 

4.4 

199.1 

57.09 

M076 

2.9 

23.0 

61.31 

M235 

20 

4.4 

199.1 

60.90 

*5““ 

1 

-1 

31.8 

•  j 

54.02 

M244 

1 

-10.1 

237.8 

53.63 

10 

-.1 

31.8 

57,45 

M244 

10 

-10.1 

237.8 

57.06 

MQ^ 

20 

-.1 

31.8 

61.26 

M244 

20 

-10.1 

237.8 

60.87 

kiflo 

Mi09 

1 

.6 

48.8 

54.18 

M255 

1 

-1.9 

268.7 

53.81 

10 

-.6 

48.8 

57.61 

M255 

10 

-1.9 

268.7 

57.24 

M109 

2Q 

-.6 

48.8 

61.42 

M255 

20 

-1.9 

268.7 

61.05 

Mm 

1 

-1.6 

^.3 

53.97 

P42A 

•• 

5.6 

71.8 

60.04 

mn 

10 

-1.6 

66.3 

57.40 

P61B 

.. 

-47.8 

264.1 

56.32 

MI38 

20 

-1.6 

66.3 

61.21 

P60A 

— 

26.3 

226.3 

60.80 

1 

3.7 

81.4 

53.88 

P59B 

— 

18.8 

88.3 

56.13 

Mioa 

10 

3.7 

81.4 

57.31 

P57A 

.. 

29.3 

127.2 

60.70 

20 

3.7 

81.4 

61.12 

P55B 

— 

29.1 

24.6 

56.51 

Table  A-2.  DISSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  6/18/90 


MLS 

Port 

DO 

Tanp 

mg/L 

“C 

DOlO 

03 

2.6 

30.5 

DOIO 

10 

3.2 

31.0 

DOlO 

20 

2.9 

29.8 

D006 

06 

4.8 

31.0 

D006 

10 

4.6 

31.9 

D006 

20 

3.7 

30.8 

M041 

01 

3.4 

31.0 

M041 

10 

3.8 

33.5 

M041 

20 

3.5 

37.0 

M066 

01 

4.9 

33.5 

M066 

10 

4.8 

33.5 

M066 

20 

4.6 

34.0 

M075 

01 

5.2 

34.0 

M075 

10 

4.9 

34.5 

M075 

20 

5.0 

36.0 

M089 

01 

4.8 

35.0 

M089 

10 

4.8 

34.0 

M089 

20 

4.3 

32.2 

M109 

01 

5.0 

32.1 

M109 

10 

4.6 

31.0 

M109 

20 

3.7 

36.0 

M138 

01 

6.0 

31.5 

M138 

10 

4.8 

32.0 

M138 

20 

5.1 

33.5 

M162 

01 

5.8 

30.5 

M162 

10 

5.0 

30.3 

M162 

20 

4.6 

29.8 

M172 

01 

5.6 

28.5 

M172 

10 

5.0 

27.8 

M172 

20 

4.8 

27.0 

MLS 

Pttrt 

DO 

Temp 

mg/L 

“C 

M181 

01 

5.4 

30.3 

M181 

10 

5.7 

28.2 

M181 

20 

5.7 

27.7 

M190 

01 

6.3 

28.0 

M190 

10 

4.9 

27.5 

M190 

20 

5.0 

27.3 

M202 

01 

6.7 

28.5 

M202 

10 

6.2 

27.7 

M202 

20 

5.0 

27.7 

M235 

01 

5.1 

31.5 

M235 

10 

5.7 

28.2 

M235 

20 

5.1 

29.5 

M244 

01 

5.4 

32.0 

M244 

10 

5.0 

31.3 

M244 

20 

5.7 

33.5 

M255 

01 

7.1 

29.0 

M255 

10 

6.1 

27.3 

M255 

20 

5.0 

27.5 

Piezmneters 


P42A  - 

3.9 

18.7 

P61B  - 

6.3 

18.9 

P60B 

2.1 

18.0 

P59B 

3.0 

19.9 

P57A  - 

5.0 

22.2 

P55B  - 

5.3 

20.0 
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T^tE  A-3.  biSSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  8/13/90 


Mis 

rutt 

PO 

mg/t 

Temp 

"C 

MLS 

Port 

DO 

mg/L 

Temp 

“C 

DOlO 

03 

5.0 

32.0 

M162 

01 

4.0 

27.0 

DOlO 

10 

5.1 

31.0 

M162 

10 

4.2 

27.0 

DOlO 

20 

3.8 

31.0 

M162 

20 

3.0 

27.0 

tiodo 

06 

4.6 

31.0 

M172 

01 

3.6 

27.0 

0006 

10 

4.5 

31.0 

M172 

10 

3.4 

27.0 

bOQ6 

20 

4.8 

29.0 

M172 

20 

2.6 

26.0 

Ni041 

01 

5.1 

32.0 

M181 

01 

4.0 

26.0 

M041 

to 

6.2 

32.0 

M181 

10 

3.2 

26.0 

M041 

20 

5.2 

31.0 

M181 

20 

3.2 

25.0 

M066 

01 

4.6 

32.0 

M190 

01 

3.8 

27.0 

M066 

id 

4.4 

31.0 

M190 

10 

3.6 

27.0 

M066 

20 

4.4 

31.0 

M190 

20 

3.4 

26.5 

M075 

01 

4.8 

30.0 

M202 

01 

5.0 

27.0 

M075 

10 

4.2 

29.0 

M202 

10 

4.0 

26.0 

M075 

20 

3.8 

28.0 

M202 

20 

2.0 

25.0 

Iild89 

01 

5.2 

28.0 

M235 

01 

5.0 

25.0 

10 

4.4 

29.0 

M235 

10 

4.6 

25.0 

M089 

3.4 

29.0 

M235 

20 

4.0 

25.0 

oi 

4-2 

28.0 

M244 

01 

4.0 

26.0 

MlO^ 

10 

3.4 

29.0 

M244 

10 

3.8 

25.0 

M109 

20 

4.6 

29.0 

M244 

20 

3.8 

26.0 

M138 

01 

4.2 

28.0 

M255 

01 

5.4 

29.0 

M13a 

10 

4.0 

28.0 

M255 

10 

5.2 

31.0 

Ml3d 

20 

3.0 

26.0 

M255 

20 

4.4 

27.0 
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TABLE  A-4.  DISSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  10/15/90 


MLS 

Port 

DO 

Temp 

MLS 

Port 

DO 

Temp 

mg/L 

“C 

mg/L 

°C 

DOlO 

03 

4.0 

31.0 

M181 

01 

4.8 

31.0 

DOlO 

10 

3.8 

30.0 

M181 

10 

4.8 

29.0 

DOlO 

15 

3.6 

31.0 

M181 

20 

5.0 

29.0 

D006 

03 

3.2 

27.0 

M190 

01 

4.6 

31.0 

D006 

10 

3.8 

28.0 

M190 

10 

4.4 

31.0 

D006 

20 

4.2 

27.0 

M190 

16 

4.2 

31.0 

M041 

01 

5.0 

30.0 

M202 

01 

5.4 

29.0 

M041 

10 

5.0 

29.0 

M202 

10 

5.2 

31.0 

M041 

15 

4.8 

30.0 

M202 

20 

4.8 

30.0 

M066 

01 

4.6 

28.0 

M235 

01 

5.2 

31.0 

M066 

10 

4.8 

31.0 

M235 

10 

4.8 

29.0 

M066 

15 

4.6 

29.0 

M235 

20 

5.0 

31.0 

M075 

01 

4.8 

29.0 

M244 

01 

5.2 

31.0 

M075 

10 

4.8 

31.0 

M244 

10 

5.0 

29.0 

M075 

16 

4.6 

28.0 

M244 

20 

4.8 

30.0 

M089 

01 

5.2 

31.0 

M255 

01 

5.6 

31.0 

M089 

10 

4.8 

30.0 

M255 

10 

5.4 

29.0 

M089 

20 

4.4 

29.0 

M255 

20 

5.6 

28.0 

M109 

01 

5.2 

31.0 

M109 

10 

4.8 

31.0 

M109 

17 

4.6 

30.0 

M138 

01 

5.2 

31.0 

M138 

10 

50 

31.0 

Piezom^ers 

M138 

20 

4.4 

31.0 

M162 

01 

4.8 

31.0 

P42A 

.. 

4.4 

20.8 

M162 

10 

4.6 

31.0 

P61B 

3.7 

19.6 

M162 

20 

4.4 

31.0 

P60A 

— 

4.6 

20.7 

M172 

01 

4.8 

31.0 

P59B 

2.6 

27.2 

M172 

10 

5.0 

31.0 

P57A 

— 

4.6 

22.2 

M172 

20 

4.8 

31.0 

P55B 

— 

5.9 

18.5 
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TABLE  A-5.  DISSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  12/4/90 


MLS 

Port 

DO 

mg/L 

Temp 

•c 

DOlO 

03 

4.0 

11.5 

nolo 

n 

3.9 

11.0 

DOlO 

15 

4.3 

10.5 

DQQ6 

08 

4.1 

12.0 

D006 

10 

4.1 

10.0 

D006 

20 

3.5 

11.0 

M(M1 

01 

4.6 

10.0 

hmi 

to 

5.0 

10.0 

M041 

19 

4.5 

12.0 

14066 

01 

5.8 

10.0 

Wi66 

10 

6.0 

11.0 

M066 

15 

4.4 

11.0 

MOTS 

01 

6.4 

12.0 

MOTS 

10 

6.T 

12.2 

MOTS 

19 

5.2 

12.5 

M009 

01 

T.O 

11.0 

MQjM 

10 

T.2 

12.0 

M0I9 

19 

7.2 

11.5 

M109 

01 

7.5 

13.5 

Hm 

10 

7.8 

14.0 

M109 

19 

6.5 

13.8 

M13S 

01 

5.0 

13.0 

M138 

10 

4.4 

13.8 

Miss 

20 

5.3 

14.5 

M162 

01 

7.2 

13.0 

Mlfi2 

10 

5.5 

14.0 

llll62 

70 

6.8 

13.5 

Mlt2 

01 

6.8 

13.0 

M1T2 

5.2 

■nm 

MLS 

Port 

DO 

Temp 

mg/L 

"C 

M181 

01 

5.5 

13.2 

M181 

10 

6.8 

14.5 

M181 

20 

6.0 

14.8 

M190 

01 

7.2 

14.5 

M190 

10 

6.6 

14.0 

M190 

20 

6.4 

14.0 

M2Q2 

01 

7.2 

16.5 

M202 

10 

7.6 

17.0 

M202 

20 

5.2 

17.5 

M23S 

01 

5.8 

16.0 

M235 

10 

7.6 

16.2 

M235 

20 

7.4 

16.0 

M244 

01 

9.0 

17.0 

M244 

10 

6.6 

16.8 

M244 

20 

6.6 

16.8 

M255 

02 

7.6 

16.2 

M255 

10 

7.7 

15.0 

M255 

20 

6.0 

14.8 

Fiesoraeters 

P42A 

3.5 

19.1 

P61B 

6.2 

17.6 

P60A 

4.3 

18.0 

P59B 

3.2 

18.2 

P<7A 

A  A 

1*7  O 

TABLE  A-6.  DISSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  3/7/91 


MLS 

Port 

DO 

Temp 

MLS 

Port 

DO 

Temp 

mg/L 

“C 

mg/L 

°C 

DOlO 

03 

2.8 

13.4 

M181 

01 

3.3 

16.3 

DOlO 

11 

2.6 

13.4 

M181 

10 

4.8 

16.1 

DOlO 

20 

3.1 

13.4 

M181 

20 

4.6 

15.6 

D006 

03 

3.5 

13.4 

M190 

01 

4.4 

16.1 

D006 

10 

2.9 

13.5 

M190 

10 

3.3 

16.0 

D006 

20 

2.6 

13.8 

M190 

20 

4.5 

16.0 

M041 

01 

3.6 

13.1 

M202 

01 

5.6 

16.0 

M041 

10 

4.0 

12.7 

M202 

10 

5.1 

15.9 

M041 

20 

3.2 

13.5 

M202 

20 

3.0 

15.7 

M066 

01 

4.5 

13.5 

M235 

01 

5.3 

16.4 

M066 

10 

4.6 

13.1 

M235 

10 

5.1 

16.1 

M066 

20 

5.2 

13.0 

M235 

20 

3.6 

16.5 

M075 

01 

5.2 

13.6 

M244 

01 

2.9 

17.1 

M075 

10 

5.3 

13.4 

M244 

10 

3.0 

16.7 

M075 

20 

7.0 

13.5 

M244 

20 

3.6 

16.5 

M089 

01 

5.2 

14.9 

M255 

01 

6.1 

16.8 

M089 

10 

5.1 

15.0 

M255 

10 

5.7 

16.5 

M089 

20 

4.3 

14.9 

M255 

20 

3.5 

16.6 

M109 

01 

4.2 

14.9 

M109 

10 

3.3 

14.9 

M109 

20 

3.5 

14.7 

M138 

01 

4.9 

14.7 

M138 

10 

3.1 

14.6 

Piezmneters 

M138 

20 

7.8 

14.3 

M162 

01 

5.0 

14.9 

P42A 

— 

7.4 

16.4 

M162 

10 

4.0 

15.0 

P61B 

— 

7.1 

17.7 

M162 

20 

4.9 

14.7 

P60A 

— 

5.0 

16.4 

M172 

01 

3.9 

14.6 

P59B 

— 

4.8 

17.3 

M172 

10 

5.0 

15.0 

P57A 

— 

8.5 

16.2 

M172 

20 

3.7 

14.6 

P55B 

6.9 

17.6 
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TABL£  A-7.  DISSOLVED  OXYGEN  AND  TEMPERATURE  DATA  ON  5/22/91 


MLS 

Jhart 

DO 

Temp 

MLS 

** - * 

rOit 

DO 

Temp 

1 

mg/L 

OC 

mg/L 

"C 

DOlO 

03 

3.8 

20.8 

M181 

01 

7.0 

26.5 

DOlO 

11 

5.2 

20.6 

M181 

10 

7.1 

26.0 

DOlO 

20 

5.0 

21.0 

M181 

20 

6.5 

26.4 

D006 

03 

4.8 

22.3 

M190 

01 

5.7 

25.0 

0006 

10 

4.7 

21.7 

M190 

10 

4.5 

24.7 

D006 

20 

4.2 

22.0 

M190 

20 

5.1 

24.8 

M041 

01 

5.0 

21.7 

M202 

03 

6.2 

24.0 

M041 

10 

5.1 

22.0 

M202 

10 

6.6 

25.7 

M041 

20 

5.0 

21.6 

M202 

20 

5.8 

25.5 

M066 

01 

6.8 

22.7 

M235 

01 

4.4 

24.3 

M066 

10 

5.2 

22.3 

M235 

10 

5.6 

23.6 

M066 

20 

5.4 

22.7 

M235 

20 

5.6 

24.2 

M075 

01 

7.2 

22.7 

M244 

01 

4.0 

24.0 

M07S 

10 

6.3 

22.3 

M244 

10 

3.5 

23.8 

i/m 

MO^ 

20 

01 

6.6 

7.8 

22.5 

24.8 

M244 

M255 

20 

01 

4.0 

7.4 

24.8 

24.2 

M089 

10 

6.4 

23.7 

M255 

10 

7.0 

24.0 

M089 

20 

5.6 

23.5 

M255 

20 

$.1 

24.6 

Mie9 

01 

5.4 

23.2 

,  M109 

10 

5.0 

23.8 

M109 

20 

4.5 

23.8 

MI3S 

01 

5.2 

23.5 

M13d 

10 

4.0 

23.2 

Piezometers 

M138 

20 

6.1 

23.0 

M1S2 

01 

4.7 

23.3 

P42A 

— 

9.2 

17.5 

MiS2 

10 

4.2 

23.3 

P61B 

— 

7.2 

18.4 

M162 

20 

4.2 

23.0 

P60A 

“ 

6.7 

18.3 

M172 

01 

6.1 

24.5 

P59B 

— 

4.3 

18.4 

M172 

10 

5.8 

24.5 

P57A 

— 

8.3 

18.3 

Uin 

20 

3.8 

23.8 

P55B 

~ 

6.5 

18.5 

i 
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SECTION  XI 
APPENDIX  B 


EVALUATION  OF  THE  SPATIAL  MOMENTS  CODE 

The  accuracy  of  the  code  used  to  estimate  the  spatial  moments  of  the  solute  plumes 
was  evaluated  through  a  set  of  numerical  tests.  In  these  tests,  three-dimensional  Gaussian 
solute  plumes  of  known  mass  were  analytically  generated  within  the  sampling  network  used 
in  the  natural-gradient  experiment.  Calculated  solute  concentrations  were  obtained  for  each 
sampling  point  in  the  network  using  Yeh’s  (1983)  analytical  solution  to  the  advection- 
dispersion  transport  equation  for  an  instantaneous  finite-volume  source.  Spatial  moments  for 
the  analytically  generated  plumes  were  then  estimated  by  the  moments  code  and  compared 
with  the  known  theoretical  values.  The  level  of  agreement  between  the  calculated  and 
theoretical  values  is  a  measure  of  the  effectiveness  of  the  particular  moments  code  that  is 
used. 

In  these  numerical  tests,  the  strength  and  dimensions  of  the  finite  volume  source  were 
similar  to  those  of  the  tritium  solution  inj«;ted  during  the  natural-gradient  study.  The  total 
injected  activity  was  0.6  Ci,  and  the  volume  of  the  source  solution  was  10  cubic  meters.  For 
an  assumed  aquifer  porosity  of  30  percent,  this  corresponded  to  a  bulk  aquifer  source  volume 
of  33.3  cubic  meters.  The  source  was  centered  at  coordinates  (0,  0,  S8.S),  and  source 
dimensions  of  4.08,  4.08,  and  2  meters  were  assigned  in  the  x,  y,  and  z  dimensions, 
respectively.  A  seepage  velocity  of  0.140  m/d  in  the  y  direction  was  assumed. 

Dispersivities  of  0.30,  0.03,  and  0.0003  meters  were  assumed  in  the  longitudinal,  lateral, 
and  vertical  directions.  These  dispersivities  produced  plumes  that  were  roughly  similar  to 
the  tritium  plumes  observed  duripg  MADE-2. 

Solute  concentration  distributions  were  computed  at  sample  points  at  elapsed  times  of 
37,  71,  142,  214,  285,  and  356  days.  Loi^tudinal  proflles  of  the  solute  concentration  for 
the  artificial  snapshots  are  shown  in  Figure  B-1.  The  overall  dimensions  of  the  artificial 
plumes  are  similar  to  those  of  the  actual  tracer  plumes  during  the  early  stages  of  MADE-2. 
However,  the  magnitude  and  variability  of  concentration  gradients  observed  in  the  field  study 
are  ncM  represented  in  the  analytically  generated  plumes. 

The  predicted  relative  activity  recovery  (M/M^)  and  plume  centroid  coordinates  (X^, 
Yc,  and  ZJ  are  given  in  Table  B-1 .  The  activity  recovery  error  oscillates  somewhat  varying 
from  -3  percent  of  the  injected  activity  at  a  displacement  of  5  meters,  to  a  maximiim  of  -1-13 
percent  at  20  meters,  and  reaches  a  minimum  of  -1  percent  at  50  meters.  These  results 
indicate  acceptable  accuracy  of  the  zero"*  moments  for  near-source  conditions  where  large 
concentration  gradients  typically  produce  the  greatest  numerical  error.  The  predicted  first 
moments  for  the  analytically  generated  plumes  are  also  reasonably  accurate.  The  estimated 
coordinates  of  the  plunw  centroid  are  within  3  percent  of  the  theoretical  mean  displacement 
for  each  case. 

These  numerical  tests  do  not  provide  a  rigorous  mathematical  verification  of  the 
mmnents  code  used.  Nevertheless,  die  results  indicate  the  code  yields  reasonably  accurate 
lower  (zero*^  and  first)  moment  estimates  for  Gaussian  plumes  defined  by  spatial 
meastttmnents  that  correspond  to  the  actual  field  sampling  network. 
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t  -  37  days 


t  »  142  days 


t  «  356  days 


OtSTANCE  FROM  INJECTXM  (2.5  X  m) 

.  LoqgibMSiisl  Profiles  of  Andytically  Oenmted  Solute  Pluiiie. 
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TABLE  B-1.  COMPARISON  OF  THE  THEORETICAL  AND  PREDICTED 

PLUME  MASS  AND  CENTROID 


Elapsed 

Plume 

Tiine(days) 

Characteristic 

Theoretical 

Predicted 

Difference 

37 

M/Mo 

1.00 

0.97 

-0.03 

Xc(m) 

0.00 

-0.04 

-0.04 

yc(ni) 

5.00 

5.15 

0.15 

Zc(m) 

58.50 

58.80 

0.30 

71 

M/Mo 

1.00 

1.11 

0.11 

Xc(m) 

0.00 

0.10 

0.10 

yc(m) 

10.00 

10.13 

0.13 

Zc(m) 

58.50 

58.50 

0.00 

142 

M/Mo 

1.00 

1.13 

0.13 

Xc(m) 

0.00 

-0.16 

-0.16 

yc(m) 

20.00 

19.89 

-0.11 

Zc(m) 

58.50 

58.50 

0.00 

214 

M/Mo 

1.00 

0.96 

-0.04 

Xc(m) 

0.00 

0.17 

0.17 

yc(m) 

30.00 

29.64 

-0.36 

Zc(m) 

58.50 

58.49 

-0.01 

285 

M/Mo 

1.00 

0.98 

-0.02 

Xc(m) 

0.00 

0.77 

0.77 

yc(m) 

40.00 

40.02 

0.02 

Zc(m) 

58.50 

58.50 

0.00 

356 

M/Mo 

1.00 

0.99 

-0.01 

Xc(m) 

0.00 

-0.04 

-0.04 

yc(m) 

50.00 

50.08 

0.08 

Zc(in) 

58.50 

58.50 

0.00 

